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ABSTRACT 


The  ability  to  use  radar  to  discriminate  ice  types  has  been  inves¬ 
tigated.  Radar  backscatter  measurements  were  made  of  shorefast  sea  ice 
near  Point  Barrow,  Alaska  in  May  1977  and  April  1978,  with  a  surface- 
based  FM-CW  scatterometer  that  swept  from  1-2  GHz  and  from  8.5-17.5  GHz. 
The  1-2  GHz  measurements  showed  that  thick  first-year  and  multiyear  sea 
ice  cannot  be  distinguished  at  10°-70°  incidence  angles,  but  that  unde- 
formed  sea  ice*can  be  discriminated  from  pressure-ridged  thick  first-year 
sea  ice’ and  lake  ice.  Results  also  indicate  that  frequencies  between 
8-18  GHz  have  the  ability  to  discriminate  between  thick  first-year  sea 
ice,  multiyear  sea  ice,  and  lake  ice.  The  lowest  frequency,  9  GHz,  was 
found  to  provide  the  greatest  separation  between  these  ice  categories 

with  significant  levels  of  separation  existing  between  angles  of  incidence 

_  D./r  (j 

from  1S9  to  7oC  The  radar  cross-sections  for  the  like  polarizations, 

VV  and  HH,  were  very  similar  in  absolute  level  and  angular  response.  The 
radar  cross-sections  for  VV-polarization  were  usually  the  highest  in 
absolute  level.  Cross-polarization  provided  a  slightly  greater  separation 
between  these  categories  of  ice.  The  8.5-17.5  GHz  angular  responses 
were  regressed  with  a  straight  l^ne\and  had  an  average  correlation  co¬ 


efficient  of  .97  (for  1%  confidence  the  correlation  coefficient  for  our 
measurements  would  need  to  be  greater  than  .89)  which  shows  the  goodness 
of  the  fit.  An  empirical  model  was  also  developed  which  describes  the 
angular  and  frequency  variations  for  thick  first-year  sea  ice,  multiyear 
sea  ice,  and  lake  ice.  The  standard  deviation  of  the  model  and  the 
averaged  scattering  cross-section  was  at  worst  2  dB.  The  individual 
measurements  had  a  worst  standard  deviation  of  3  dB  about  the  average 
scattering  cross-section. 
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1.0  INTRODUCTION 


The  Arctic  is  the  area  of  the  far  north  characterized  by  distinc¬ 
tively  polar  conditions  of  climate,  plant  life,  etc.  The  term  Arctic 
is  derived  from  the  Greek  word  "arktos"  meaning  "bear"  which  refers  to 
the  northern  constellation  of  the  Bear.  International  interest  in  the 
Arctic  and  sub-Arctic  regions  has  steadily  increased  during  the  twen¬ 
tieth  century  and  particularly  since  the  Second  World  War.  Three 
major  factors  have  led  to  this  rising  interest:  the  advantage  of 
the  North  Pole  route  as  a  short-cut  between  important  centers  of  popu¬ 
lation;  the  growing  realization  of  economic  potentialities  such  as 
minerals  (especially  petroleum),  forest  resources,  and  grazing  areas; 
and  the  importance  of  the  regions  in  the  study  of  global  meteorology. 

The  Arctic  Ocean  is  partially  covered  with  sea  ice  at  all  times. 
With  a  surface  area  estimated  at  5,427  square  miles,  ice  may  have  an 
extent  ranging  from  1/10  cover  to  9/10  or  10/10  cover,  depending  upon 
the  season  and  locality.  The  ice  cover,  pack  ice,  is  not  a  static 
homogeneous  film  upon  the  sea  surface.  It  consists  of  large  and  small 
areas  of  open  water,  areas  with  pressure  ridges,  rafted  ice,  etc., 
forming,  growing,  and  melting  and  consequently  covers  areas  of  variable 
extension. 

Due  to  the  increasing  variety  of  scientific  and  operational 
interests  in  the  Arctic  polar  region,  testing  and  evaluation  of  state-' 
of-the-art  remote  sensors  has  been  in  progress  for  many  years  to 
determine  their  potential  capability  in  sea  ice  studies.  Development 
of  microwave  frequency  radars  used  as  remote  sensors  broadens  pros¬ 
pects  of  remote  measurements  of  such  interesting  and  dynamic  features 
of  the  Earth's  surface.  The  capability  of  radar  to  provide  useful 
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information  under  weather  and  lighting  conditions  unsuitable  for  visual 
sensing  attracts  considerable  interests  toward  developing  it  into  a 
tool  to  be  used  in  sea  ice  reconnaissance  and  forecasting. 

1.1  Background 

Study  and  exploration  of  the  Arctic  probably  began  with  a  Greek 
named  Pytheas,  who  in  the  fourth  century  BC  made  an  astonishing  voyage 
from  the  Mediterranean  around  Britain  and  reached  a  place  he  called 
Thule,  variously  identified  as  the  Shetlands,  Iceland,  and  Norway. 

His  accounts  of  this  remarkable  voyage  were  for  centuries  discredited, 
but  the  idea  of  his  Thule,  shrouded  in  fog  and  believed  to  be  the  end 
of  the  Earth,  caught  the  imagination  of  many.  From  these  first 
glimpses  of  the  inhabitable  "frigid  zone"  shrouded  in  a  mist  of  super¬ 
stitious  beliefs,  primitive  navigation  techniques,  and  an  obscurity 
resulting  from  inaccurate  ideas  of  the  shape  of  the  Earth,  exploration 
continued  by  men  who  could  not  quench  their  thirst  for  adventure  by 
more  calm  and  rudimentary  methods. 

A  comprehensive  expansion  of  all  types  of  research  in  the  Arctic 
Ocean  took  place  between  1950  and  1970.  Many  drift  ice  research  sta¬ 
tions  were  established  by  the  United  States  and  the  Soviet  Union 
serving  as  bases  for  studies  of  atmosphere,  ice,  ocean  waters,  and 
crust  beneath.  U.S.  nuclear  submarines  made  their  first  voyages 
through  the  Arctic  Ocean  in  1957,  making  possible  continuous  profiles 
of  many  geophysical  parameters.  Research  efforts  by  the  U.S.  recently 
culminated  in  the  Arctic  Ice  Dynamics  Joint  Experiment  (a  study  of 
basic  processes  and  data  collecting  for  use  in  the  modeling  of  sea  ice). 


The  naval  military  importance  of  the  Arctic  has  risen  significantly 
over  the  past  few  years.  The  chain  of  events  leading  to  policy  of 
energy  independence  and  the  consequent  emphasis  on  accelerated  explora¬ 
tion  of  our  North  Slope  and  offshore  Alaskan  oil  reserves  has  served  as 
one  cause.  As  a  consequence,  strategists  have  considered  the  vital 
role  of  our  naval  forces  in  the  protection  of  a  new  and  important  sea 
lane.  A  second  factor  has  been  the  increasing  long-range  capabilities 
of  submarine-launched  ballistic  missiles,  making  the  Arctic  Ocean  a 
possible  patrol  and  launch  area.  A  third  factor  involves  freedom  of 
the  seas  and  geographical  intentions.  Very  recently  countries  bordering 
the  Arctic  have  threatened  to  extend  their  rights  over  contiguous 
waters  to  the  200  mile  limit,  which  could,  in  some  measures,  bottle- 
up  the  narrow  eastern  entrance  to  the  Arctic  Basin.  The  nuclear  attack 
submarine,  with  its  exceptional  mobility  in  ice-infested  waters,  has 
an  important  potential  role  in  all  of  the  above  considerations. 

The  polar  regions  greatly  effect  the  environment  of  the  Earth's 
surface.  In  the  global  thermodynamics  cycle  of  atmosphere  and  ocean, 
the  polar  regions  are  the  heat  sinks.  In  the  period  of  one  revolution 
of  the  Earth,  the  polar  regions  lose  more  heat  to  space  than  is  received 
from  the  sun.  As  a  result,  they  are  regions  which  are  cold  and  maintain 
permanent  ice  cover  of  annually  ranging  extent.  In  compensation  for 
the  heat  lost  to  space,  heat  from  the  lower  latitudes  is  circulated  into 
the  polar  regions.  The  vertical  extent  of  sea  ice  cover  is  extremely 
small  compared  to  its  vast  horizontal  extent  (about  1:10^).  Small 
perturbations  in  the  heat  balance  in  the  ice  surface  may  result  in 
large  alterations  in  the  ice  cover,  resulting  in  changes  of  terrestrial 


4. 

albedo,  sea  surface  temperature,  ocean  mixing,  evaporation,  and  so 
forth. 

Sea  ice  is  a  dynamic  obstacle  of  shipping.  For  global  weather 
forecasting  it  is  a  prized  piece  of  the  puzzle,  an  integral  part  of 
the  global  weather  model.  On  a  large  time  scale  ice  is  forming, 
growing,  and  melting;  consequently  it  covers  areas  of  variable  exten¬ 
sion.  On  a  short  time  scale  weather,  mainly  the  wind,  influences  the 
ice  cover  causing  movement  of  the  ice,  formation  of  open  leads  and 
ice  ridges,  etc.  In  order  to  make  navigation  as  safe  as  possible  and 
to  optimize  it  economically  by  choosing  the  easiest  routes  through 
the  ice,  it  is  necessary  to  watch  the  ice  situation.  Ice  cover,  being 
a  key  in  global  weather,  must  be  monitored.  This  means  mapping, 
analyzing,  and  forecasting  the  ice  situation. 

For  sea  ice  mapping  may  be  the  name  of  the  game.  Large  areas  of 
sea  ice  have  been  mapped  in  the  past,  primarily  by  the  use  of  conven¬ 
tional  aerial  photographic  methods  using  black  and  white,  color,  and 
infrared  films.  It  was  shown  by  Anderson  (1970)  that  ages  and  discrim¬ 
ination  among  categories  of  sea  ice  could  be  determined  from  stereo¬ 
scopic  aerial  photography  by  employing  the  usual  photo-interpretation 
techniques  of  terrain  air-photo  analysis.  Pattern  configuration,  sur¬ 
face  features,  and  tone  are  used  in  the  interpretation.  Local  weather 
and  daytime  usage  limit  the  operational  capability  of  this  remote 
sensor  in  the  uncertain  weather  conditions  and  durations  of  inadequate 
incident  light  characteristic  of  the  Arctic  region. 

Development  of  remote-sensing  radar  systems  operating  in  the  micro- 
wave  region  of  the  frequency  spectrum  has  broadened  the  prospects  of 
remotely  measuring  significant  features  of  the  Earth's  surface.  Since 


the  performance  of  the  remote-sensing  radars  is  unaffected  by  poor 
weather  conditions  and  lack  of  incident  light,  the  application  of 
radar  in  the  remote  sensing  of  sea  ice  appears  to  be  particularly 
attractive. 

Tie  potential  application  of  SLAR  to  sea  ice  research  and  recon¬ 
naissance  was  recognized  in  1962,  when  the  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory  utilized  a  U.S.  Air  Force  AN/APQ-56 
K-band  SLAR  in  sea  ice  experiments  conducted  during  April  and  August  of 
that  year  [1].  This  study  indicated  that  various  important  sea  ice 
features  could  be  identified  on  the  K-band  imagery,  and  that  the  SLAR 
system  provided  the  capability  of  obtaining  good-resolution  image  maps 
of  large  expanses  of  ice  independent  of  incidence-light  and  weather 
conditions. 

Since  that  time  radar  has  been  used  extensively  to  monitor  sea 
ice  by  Rouse  (1969);  Johnson  and  Farmer  (1971);  Glushkov  and  Komarov 
(1971);  Ketchum  and  Tooma  (1973);  Parashar  et  al .  (1974);  Dunbar  (1975); 
Dunbar  and  Weeks  (1976);  Gray  et  al .  (1977);  and  Ketchum  (1977). 

These  studies  have  primarily  been  made  at  2  and  3  cm  wavelengths, 
although  a  few  studies  were  made  at  25  cm  wavelengths.  Most  SLAR 
flights  have  used  angles  of  incidence  near  grazing.  These  studies 
ter  to  agree  that,  with  some  reservations,  the  ability  to  detect  and 
interpret  ice  age,  ice  drift,  ice  surface  topography  such  as  pressure- 
ridges  and  fractures,  and  to  evaluate  sea  ice  conditions,  in  general, 
was  very  good. 

The  first  attempt  to  expand  the  knowledge  of  the  parametric 
response  of  radar  systems  sensing  sea  ice  began  in  1967  when  the 
National  Aeronautics  and  Space  Administration,  Navy  Oceanographic 
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Office,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory, 
Arctic  Institute  of  North  America,  and  the  University  of  Kansas  studied 
the  ability  of  a  2.25  cm  wavelength  radar  scatterometer  to  identify 
different  ice  types.  The  radar  scatterometer  is  a  calibrated  instrument 
designed  to  measure  the  differential  scattering  coefficient  a0  as  a 
function  of  incidence  angle  [11].  Analysis  of  the  1967  scatterometer 
data  was  carried  out  by  Rouse  [2],  who  demonstrated  that  a  2.25  cm 
wavelength  (13.3  GHz)  vertically  polarized  scatterometer  could  be  used 
to  discriminate  ice  types.  For  instance,  multiyear  returns  were 
higher  than  thick  first-year  ice  at  all  angles  of  incidence. 

In  April  1970  another  joint  experiment  was  conducted  by  NASA, 

Naval  Oceanographic  Office,  and  the  University  of  Kansas  in  which 
systematic  radar  backscatter  measurements  of  sea  ice  were  made  at  400 
MHz  (HH,  VV,  VH,  and  HV  polarizations)  and  13.3  GHz  (VV  polarization). 
The  results  showed  that  multiyear  ice  gave  the  strongest  return  at 
13.3  GHz  [12],  First-year  ice  and  open  water  gave  the  strongest 
return  at  400  MHz.  Also  discrimination  of  ice  from  open  water  was 
possible  with  both  frequencies.  Four  polarization  16.5  GHz  radar 
imagery  was  also  analyzed.  Results  showed  that  0-18  cm,  18-90  cm, 

90-360  cm  ice  and  open  water  were  discriminate.  These  results  were 
consistent  with  the  scatterometer  results.  Cross-polarization  appeared 
to  be  better  in  discriminating  sea  ice  types,  but  problems  with  the 
images  prevented  a  firm  conclusion. 

During  the  winter  and  spring  of  1975-1976,  measurements  were  made 
with  a  mul ti -polarized  13.3  GHz  scatterometer  by  the  Canada  Centre 
for  Remote  Sensing  in  conjunction  with  Aidjex  and  Beaufort  Sea  offshore 
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programs  [9].  Results  showed  systematic  changes  in  microwave  back- 
scatter  which  strongly  correlated  with  gross  ice  type  categories. 
Multiyear  ice  showed  significantly  higher  backscatter  than  thick  first- 
year  ice.  Typical  differences  in  backscatter  ran  8-10  dB  for  like- 
polarization  and  15-18  dB  for  cross-polarization. 

An  X-band  SLAR  experiment  was  also  conducted  recently  (May  1975) 
by  the  Naval  Oceanographic  Office  over  areas  of  the  Lincoln  Seas, 

Nares  Strait,  and  Baffin  Bay.  The  ability  to  interpret  sea  ice 
features  and  evaluate  sea  ice  conditions  was  considered  to  be  generally 
good  in  the  more  northernly  Arctic  areas.  However,  the  young  stages 
of  ice  development  observed  in  Baffin  Bay  made  interpretation  and 
discrimination  much  more  difficult  due  to  a  high  homogeneous  radar 
return.  Multiyear  ice  was  easy  to  discriminate  due  to  predictable 
return  and  the  well-defined  shapes  of  the  floes.  Younger  ice  was 
discriminable  from  old  ice,  in  a  gross  sense,  but  discrimination  among 
the  several  stages  of  young  ice  development  was  not  possible  from 
backscatter  information  and  pattern  analysis.  New  ice-ridges  and 
hummocks  provided  strong  returns.  Old  weathered  ridges  and  hummocks 
provided  backscatter  returns  that  were  similar  to  returns  observed  in 
undeformed  old  ice  suggesting  that  volume  scattering  is  the  major 
backscatter  mechanism  of  multiyear  ice. 

1 . 2  Review  of  Previous  Results 

Previous  scatterometry  measurements  of  sea  ice  have  been  detailed 
in  papers  by  Rouse  (1968),  Parashar  et  al.  (1974),  Gray  et  al.  (1977), 
and  Gray  et  al.  (1977).  Each  of  these  will  be  considered  individually, 


7. 


with  brief  descriptions  given  of  the  radars  used,  the  sea  ice  conditions 
at  time  of  measurement,  and  their  results. 
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1.2.1  J.W,  Rouse  (1968). 

Radar  backscatter  measurements  were  made  at  13.3  GHz  (VV  polarization). 
These  measurements  were  made  using  the  NASA  Ryan  Redop  scatterometer. 

The  illuminated  area  subtended  120°  in  the  along-track  direction  and  3° 
in  the  across-track  direction.  The  measurements  were  made  off  the  coast 
of  Alaska  near  Point  Barrow  during  the  month  of  May  in  1967.  Ground 
truth  was  provided  by  black-and-white  photography  and  color  infrared 
photography.  Some  documentation  was  done  by  ground  parties. 

The  results  (of  flight  line  92)  shown  in  Figure  1.2-1  are  those 
for  an  old  multiyear  ice  floe  and  smooth  first-year  ice.  The  multiyear 
ice  floe  was  found  to  be  several  years  old  and  was  well -weathered. 
Hummocks  were  well-rounded.  The  smooth  thick  first-year  ice  had  an 
average  thickness  of  1.35  meters,  was  uniform,  and  was  covered  with  1-2 
cm  of  snow.  Multiyear  ice  returns  were  found  to  be  greater  than  5  dB 
higher  than  first-year  ice  returns  at  angles  greater  than  25°.  Figure 
1.2-2  shows  the  return  from  another  flight  line  (line  94).  The  ice 
types  and  conditions  were  found  to  be  similar  to  the  previous  line. 

This  figure  shows  the  returns  from  the  multiyear  ice,  ridged  ice,  and 
first-year  ice.  Figure  1.2-3  shows  the  returns,  obtained  from  a  third 
flight  line  (line  91),  of  smooth  first-year  ice,  ridged  first-year  ice, 
and  open  water.  These  results  showed  that  it  was  difficult  in  this 
case  to  discriminate  between  smooth  and  ridged  first-year  ice.  It 
was  also  demonstrated  that  open  water  was  discriminable  from  first-year 
ice.  Figure  1.2-4  shows  the  variation  in  the  return  of  first-year  Ice 


for  the  three  flight  lines,  with  good  agreement  in  shape.  It  was  noted 
that  the  ice  in  line  91  was  very  smooth  compared  to  the  block-strewn 
ice  of  the  other  lines.  Figure  1.2-5  shows  the  variation  in  the  multi¬ 
year  ice  returns  found  in  flight  lines  92  and  94. 

1.2.2  S.K.  Parashar,  A.W.  Biggs,  A.K.  Fung  and  R.K.  Moore  (1974). 
Radar  backscatter  measurements  of  sea  ice  were  made  at  400  MHz  (VV,  HH, 
VH,  HV  polarizations)  and  13.3  GHz  (VV  polarization).  The  NASA  Ryan 
Redop  scatterometer  was  used  in  the  13.3  GHz  measurement  program.  The 
illuminated  area  subtended  120°  in  the  along-track  direction  and  3°  in 
the  across-track  direction.  The  NASA-Emerson  scatterometer  system  was 
used  in  the  400  MHz  program.  The  illuminated  area  subtended  120°  in 
the  along-track  direction  and  7.5°  in  the  across-track  direction.  This 
mission  was  conducted  off  the  coast  of  Alaska  near  Point  Barrow  during 
the  month  of  April  1970.  The  surface  truth  was  provided  by  stereo 
aerial  photographs.  The  ice  surface  temperature  of  -10°  C  (obtained 
from  a  surface  party)  was  assigned  to  all  categories  of  ice.  Seven 
categories  were  used  to  describe  the  ice  conditions.  These  were: 


Category  1 :  Open  water 
Category  2:  New  ice  (0-5  cm  thick) 

Category  3:  Thin  young  ice  (5-18  cm  thick) 

Category  4:  Thick  young  ice  (18-30  cm  thick) 
Category  5:  Thin  first-year  ice  (30-90  cm  thick) 
Category  6:  Thick  first-year  ice  (90-180  cm  thick) 
Category  7:  Multiyear  ice  (180-360  cm  thick) 


Results  with  the  13.3  GHz  scatterometer  showed  that  open  water 
gave  the  strongest  return  at  near-vertical  angles,  whereas  beyond  15° 


INCIDENCE  ANGLE  IN  DEGREES 


Comparison  of  scattering  coefficient  variations  for  multiyear 
ice  (Rouse,  1968). 
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multiyear  ice  gave  more  return  than  the  other  categories  of  ice  or  water. 
Thick  first-year  ice  was  the  next  strongest  in  return  at  these  angles 
and  was  4  dB  less  strong  than  multiyear  ice  returns  (see  Figure  1.2-6). 

At  the  larger  angles  both  multiyear  ice  and  thick  first-year  ice  were 
discriminable  from  the  younger  stages  of  ice  and  from  open  water. 

Also,  at  these  angles  there  was  some  ability  to  discriminate  between 
new  and  thin  young  ice  and  thin  first-year  and  thick  first-year  ice. 

Radar  cross-sections  versus  thickness  curves  were  produced  to  demon¬ 
strate  the  ability  to  discriminate  among  ice  types.  These  curves  are 
shown  in  Figure  1.2-7.  As  these  curves  demonstrate,  the  ability  to 
discriminate  is  a  function  of  the  angle  of  incidence.  Also,  there  are 
ambiguities  in  the  radar  returns  from  ice  under  18  cm  thick  and  ice 
90-180  cm  thick.  This  ambiguity  exists  because  the  thinner  ice  gives 
more  return  than  the  moderately  thick  ice. 

The  results  with  the  400  MHz  scatterometer  for  VV  and  VH.  polari¬ 
zations  are  shown  in  Figures  1.2-8  and  1.2-10,  respectively.  For  VV 
polarization  open  water  gave  the  strongest  return  at  all  angles.  New 
and  thin  young  ice  gave  the  minimum  return  at  all  angles  and  thin  and 
thick  first-year  ice  gave  more  return  than  multiyear  ice.  For  HH 
polarization,  open  water  returns  were  the  strongest  at  angles  less 
than  20°,  whereas  thin  and  thick  first-year  ice  produced  the  strongest 
returns  beyond  20°.  Again,  ice  less  than  18  cm  thick  gave  minimum 
returns  at  all  angles.  For  both  cross-polarizations  (HV  and  VV),  open 
water  produced  at  almost  all  angles  the  strongest  returns.  Again,  ice 
less  than  18  cm  thick  produced  the  minimum  returns.  Radar  cross- 
section  versus  ice  thickness  curves  are  shown  in  Figures  1.2-9  and 
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Figure  1.2-8.  Experimental  a0  versus  9  for  different  categories,  400  MHz, 
VV  (Parashar  et  al . ,  1974). 
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Figure  1.2-9.  Experimental  3°  versus  ice  thickness  for  different  9,  400  MHz 
VV  (Parashar  et  al.,  1974). 


1.2-11 ,  respectively.  The  dynamic  range  in  scattering  cross-section 
was  found  to  be  greater  with  cross-polarization  than  with  like-polar¬ 
ization. 

1.2.3  L.  Gray,  J.  Cihlar,  and  S.  Parashar  (1977). 

Radar  backscatter  measurements  of  sea  ice  were  made  at  13.3  GHz  (HH 
and  HV  polarizations).  These  measurements  were  made  using  the  Ryan 
Model  720  dual-polarized  scatterometer  system.  The  illuminated  area 
subtended  120°  in  the  along-track  direction  and  3°  in  the  across-track 
direction.  The  measurements  were  made  off  the  east  coast  of  Canada 
during  the  winter  and  spring  of  1975-1976.  The  primary  ground  truth 
was  provided  by  a  RC-10  camera  and  a  Daedalus  infrared  scanner.  In- 
situ  measurements  of  temperature,  salinity,  snow  cover,  and  thickness 
of  sho refast  ice  were  made. 

An  area  called  Forteau  Bay  was  investigated.  Measurements  were 
made  of  areas  containing  smooth,  medium  rough,  and  rough  first-year 
ice,  rafted  young  ice,  new  ice,  and  open  water.  Results  showed  that, 
as  expected,  open  water  gives  the  maximum  return  near  vertical  inci¬ 
dence  and  then  decreases  relatively  rapidly  with  increasing  angle  (see 
Figures  1.2-12  and  1.2-13).  For  this  group  of  ice  types,  rough  first- 
year  ice  gives  maximum  return  at  angles  greater  than  20°,  followed  by 
medium  rough  first-year  and  then  smooth  first-year  ice.  Rafted  young 
ice  was  found  to  have  an  angular  response  very  similar  to  the  response 
of  medium  rough  thick  first-year  ice,  even  though  the  rafter  ice  was 
thinner  than  the  thick  first-year  ice.  It  is  important  to  note  that 
a  range  of  at  least  20  dB  was  noted  in  the  radar  return  from  the 
three  shorefast  ice  regions  which  were  described  by  their  surface 


roughness.  The  results  for  cross-pol arization  were  similar  to  those 
obtained  for  like-polarization,  although  the  level  of  separation 
between  ice  type  categories  appeared  to  be  slightly  greater  in  the 
case  of  cross-polarization. 

An  area  called  Northumberland  Strait  was  also  investigated. 
Measurements  were  made  of  the  following  ice  categories:  thin  first- 
year,  consolidated  brash,  rafted  young,  heavily  rafted  nil  as,  rafted 
nilas,  and  open  water.  The  angular  response  of  the  scattering  coeffi¬ 
cients  are  shown  in  Figures  1.2-14  and  1.2-15.  Among  these  ice  types 
thin  ice  gave  the  strongest  return  throughout  all  angles.  Rafted 
nilas  produces  a  relatively  low  return  in  comparison  to  other  ice 
returns,  whereas  heavily  rafted  nilas  produced  returns  that  were 
slightly  higher  than  those  of  rafted  young  ice.  The  order  of  return 
from  these  various  ice  types  for  cross-polarization  is  similar  to  like- 
polarization  returns,  but  separation  between  ice  categories  is  slightly 
greater  than  for  like-polarization.  In  this  case,  rafted  nilas 
produced  returns  that  were  significantly  lower  than  return  from  open 
water. 

1.2.4  Gray,  Ramseier,  and  Campbell  (1977). 

Radar  backscatter  measurements  were  made  at  13.3  GHz  (HH,  VV,  HV,  VH 
polarizations).  These  measurements  were  made  using  a  Ryan  Model  720 
scatterometer  system.  The  illuminated  area  subtended  120°  in  the  along- 
track  direction  and  3°  in  the  across-track  direction.  The  measurements 
were  made  in  the  Beaufort  Sea  during  the  month  of  March  in  1975. 

Ground  truth  was  provided  by  black-and-white  photography  and  infrared 
imagery.  Detailed  dielectric  constant,  temperature  and  salinity  pro¬ 
files  were  obtained  for  a  variety  of  ice  types  by  surface  parties. 
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Figure  1.2-12.  Scattering 
coefficient  a0  versus  angle  9 
Forteau  Bay,  13.3  GHz  HH  (Gray 
et  al.,  1977). 


Figure  1.2-13.  Scattering  coefficient 
o°  versus  angle  0  Forteau  Bay,  13.3 
GHz,  HV  (Gray  et  al . ,  1977). 


Figure  1.2-14.  Scattering  coefficient  a0  versus 
angle  9  Northumberland  Strait,  13.3  GHz,  HH 
(Gray  et  al . ,  1977) . 
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Figure  1.2-15.  Scattering  coefficient  a0  versus 
angle  e  Northumberland  Strait, 
13.3  GHZ,  HV  (Gray  et  al . ,  1977). 


Figure  1.2-16.  Scattering  coefficient  o°  versus 

angle  8  Bay  of  Chaleur,  13.3  GHz  HH 
(Gray  et  al . ,  1977). 


Their  results  indicated  that  there  was  a  strong  correlation  between 
the  magnitude  of  the  radar  scattering  cross-section  and  major  ice  types. 

Old  ice  gave  significantly  higher  backscatter  returns  than  first-year 
and  young  ice  forms.  Also,  at  all  angles  the  average  return  for  the 
old  ice  is  greater  than  the  maximum  signal  for  the  ridges  on  the  first- 
year  ice.  The  angular  response  for  an  old  ice  floe  is  shown  in  Figure 

1.2- 17.  Scattering  curves  for  smooth  young  ice  are  shown  in  Figure 

1.2- 18.  The  scattering  curves  for  first-year  ice  is  shown  in  Figure 

1.2- 19.  The  error  bars  that  are  plotted  on  the  curves  indicate  +  1  o 
spread  in  the  results. 

1.2.5  Comparison  of  Backscatter  Returns  at  13.3  GHz  (VV-polarization). 
The  angular  responses  of  multiyear  and  thick  first-year  ice  obtained  by 
Rouse,  Parashar  and  Gray  are  presented  in  Figure  1.2-20.  In  each  case, 
at  incidence  angles  greater  than  10°,  multiyear  ice  returns  were  higher 
than  thick  first-year  ice  returns.  Average  separations  between  respon¬ 
ses  for  thick  first-year  and  multiyear  ice  were  6.0  dB  for  Rouse,  3.5  dB 
for  Parashar,  and  12.0  dB  for  Gray.  The  angular  responses  of  multiyear 
ice  in  the  cases  of  Rouse  and  Gray  agreed,  both  in  level  and  in  shape, 
up  to  incidence  angles  of  40°.  The  angular  responses  of  thick  first- 
year  ice  in  the  cases  of  Rouse  and  Gray  were  somewhat  similar  in  level 
and  shape  at  near-vertical  angles  of  incidence.  Significant  level 
differences  (10  dB)  existed  between  the  results  of  Parashar  and  the 
results  of  Rouse  and  Gray.  Table  1.2-1  provides  for  a  comparison  of 
radar  cross-section  levels. 


Absolute  levels  and  the  angular  response  of  experimentally  measured 
radar  scattering  cross-sections  of  a  target  may  vary  from  sensor  to  sen¬ 
sor  as  was  demonstrated  in  the  above  results.  These  variations  tend  to 
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2-17.  Graph  illustrating  the  four  backscattering  curves  as  a 
function  of  incidence  angle  for  relatively  smooth, 
ridge-free  old  ice  (Gray,  Ramseier,  and  Campbell, 
1977). 


KACKSCATTER^G  COEFFICIENT 
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scattering  curves  for  smooth  young  ice 

AIOOEX,  8EAUFORT  SEA,  MARCH  1975 


Figure  1. 


2-18.  Graph  illustrating  the  four  backscattering  curves  as 
function  of  incidence  angle  for  an  area  of  smooth 
young  ice  (Gray,  Ramseier,  and  Campbell,  1977). 


SCATTERING  CURVES  FOR  FIRST  TEAR  ICE 
6I0JEX,  BEAUFORT  SEA,  MARCH  1975 


Figure  1.2-19.  Graph  illustrating  the  four  backscattering  curves  as  a 
function  of  incidence  angle  for  an  area  of  first-year 
ice  (Gray,  Ramseier,  and  Campbell,  1977). 


Radar  Cross  Section  -o°  (dB) 


Angle  of  Incidence  (Degrees) 


Figure  1.2-20.  Comparison  of  multiyear  and  thick  first-year  ice 

scattering  cross-sections  according  to  Rouse  (1968), 
Parashar,  et  al .  (1974),  and  Gray,  et  al .  (1977). 
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arise  due  to  two  problem  areas  that  are  common  to  airborne  scatterometer 
systems.  First,  there  is  the  problem  of  accurately  describing  the  illum¬ 
ination  on  the  ground.  To  do  this  properly  a  gain  pattern  of  the 
antenna  must  be  measured  with  the  antenna  attached  to  a  very  good  rep¬ 
lica  of  the  fuselage  of  the  airborne  platform.  Only  in  this  way  is  it 
really  known  for  certain  how  the  gain  distribution  characteristics  are 
shaped  by  the  multi  curved  shape  of  the  airborne-ground  plane.  Second 
unless  extreme  care  and  precision  is  taken  with  the  calibration  of 
the  airborne  sensor  with  a  known  standard  radar  target,  and  is  done 
periodically,  variations  in  the  absolute  level  in  calibration  may  occur. 


TABLE  1.2-1 

COMPARISON  OF  THE  RADAR  CROSS-SECTION  o°  (IN  dB)  RESULTS 
ROUSE  (1968),  PARASHAR  (1974),  AND  GRAY  (1977)  FOR 
MULTIYEAR  AND  THICK  FIRST-YEAR  SEA  ICE  RADAR  CROSS-SECTIONS 


MULTIYEAR 


THICK  FIRST-YEAR 


Rouse 

Parahsar 

Gray 


o 

o 

o 

o 

CM 

60° 

20° 

o 

O  I 

1 

60° 

-7.0  -13.5 

-19.0 

-13.0 

-19.0 

-26.0 

+  .5  -  1.0 

-12.5 

-  3.5 

-  4.5 

-  6.0 

-7.0  -11.5 

-13.5 

-15.0 

-24.0 

-27.5 

2.0  CHARACTERISTICS  OF  SEA  ICE 


2.1  Formation 

Sea  ice  forms  due  to  freezing  of  sea  water.  Sea  ice  consists  of 
fresh  water  ice  and  very  small  liquid  inclusions  of  brine,  referred  to 
as  brine  pockets,  entrapped  in  the  ice  matrix  during  the  process  of 
freezing.  Fresh  water  ice  forms  at  temperatures  below  the  initial 
freezing  point  temperature,  at  which  ice  is  in  equilibrium  with  the 
sea  water  of  given  salinity.  The  impurities  present  in  the  sea  water 
are  rejected  during  the  formation  of  pure  ice,  enriching  the  neighboring 
liquids.  This  is  why  the  salinity  of  sea  ice  at  all  stages  of  develop¬ 
ment  is  less  than  the  original  sea  water  from  which  it  was  formed. 

Many  factors  influence  the  formation  of  sea  ice.  These  include 
wind,  sea  state,  currents,  cooling  rate,  distribution,  and  salinity 
of  the  sea  water,  surface  temperature,  and  depth  of  water.  As  has 
been  found,  the  combination  of  these  factors  lead  to  variation  in  the 
vertical  and  horizontal  properties  throughout  the  ice  sheet.  Also, 
under  conditions  of  wind,  current,  waves,  and  pressure,  ice  fields 
may  break  up  into  floes  and  drift.  This  redistribution  may  cause 
pressure  effects  which  may  cause  ice  to  pile  up  as  ridges  and  hummocks 
or  the  overriding  of  one  sheet  of  ice  over  another,  which  is  called 
rafting. 

Sea  ice  may  be  classified  according  to  its  age,  the  way  it  was 
formed,  or  by  concentration.  A  widely  used  system  of  classification 
is  the  World  Meteorological  Organization  "Glossary  of  Ice  Terms"  [19]. 
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2.2  Physical  Properties 

The  ability  to  distinguish  between  open  water,  thin  ice,  first- 
year  ice,  and  multiyear  ice  is  of  prime  importance.  For  example, 
thickness  and  physical  properties  differ  greatly  between  first-year 
and  multiyear  ice.  In  general,  a  first-year  ice  floe  is  smooth-sur¬ 
faced  with  angular  floe  edges.  Multiyear  ice  has  generally  a  rolling 
rough  surface  while  the  edges  of  the  floe  are  rounded.  These  changes 
occur  due  to  the  melt  conditions  that  take  place  in  the  summer  which 
modify  surface  features  and  the  grinding  action  that  takes  place 
between  floes  tend  to  round  off  the  corners.  The  melt  process  not 
only  modifies  the  surface  features,  but  extensively  changes  the  material 
properties  of  the  ice.  Typical  examples  of  the  salinity,  density,  and 
temperature  profiles  of  first-year  and  multiyear  ice  are  illustrated 
in  Figure  2.2-1.  In  first-year  ice  the  salinity  is  usually  high, 
both  near  the  air-ice  interface  and  at  the  ice-water  interface.  The 
density  may  be  slightly  lower  near  the  air-ice  interface,  but  is  con¬ 
stant  throughout  the  remainder  of  the  ice  sheet.  The  temperature  pro¬ 
file  is  very  often  almost  linear  during  most  of  the  year.  If  the  ice 
is  not  covered  by  a  layer  of  snow,  the  surface  temperature  of  the  ice 
is  close  to  the  air  temperature.  Because  the  thermal  conductivity  of 
snow  is  much  lower  than  that  of  ice,  a  layer  of  snow  upon  the  ice 
keeps  the  surface  temperature  of  ice  closer  to  the  melting  point. 

Bare  sea  ice  temperatures  track  fairly  well  with  rapidly  changing 
air  temperatures .  Snow  cover  tends  to  keep  the  temperature  of  the 
ice  constant  with  changes  of  10°  C  in  air  temperature  needed  to 
change  the  temperature  of  the  ice  surface  by  5°  C  with  these  tempera¬ 
ture  fluctuations  effecting  the  top  25  cms  or  so.  In  multiyear  ice 


r 


2.2-1-  Depth  profiles  of  salinity,  density,  and  temperature  for 
first-year  and  multiyear  ice  obtained  during  April  1975 
at  the  AIDJEX  camo.  "he  schematic  of  the  core  is  shown 
on  the  left  hand  side  of  the  orofiles  indicating  the 
structural  changes  in  the  ice  cover  (Camobell  et  a!.,  IS 
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the  salinity  starts  near  zero  and  increases  with  depth.  Density  is 
lower  in  the  top  part  of  the  cover.  The  temperature  profile  follows  in 
the  same  fashion  as  in  first-year  ice. 

The  phase  diagram  for  standard  sea  ice  may  be  used  to  describe 
some  of  the  physical  properties  of  sea  ice  (see  Figure  2.2-2).  The 
amount  of  brine  present  may  be  found  from  the  diagram.  As  the  tempera¬ 
ture  is  lowered,  the  brine  concentration  increases  as  the  volume  of 
brine  decreases.  As  shown  by  the  phase  diagram,  sea  ice  at  equilibrium 
at  a  given  temperature  consists  of  fresh  water  ice  and  brine  of  a 
specific  composition.  The  first  salt  to  precipitate  is  NaSO^  •  10  HjO, 
which  begins  to  crystallize  at  -8.2°  C.  Brine  concentration  continues 
until  -22.9°  C  when  NaCl  •  2  HjO  starts  to  precipitate.  NaCl  ,  being 
the  principal  salt  in  the  brine  liquid,  causes  a  rapid  decrease  in 
brine  volume  near  -23°  C.  8rine  liquid  exists  at  temperatures  below 
-50°  C. 

The  relative  volume  of  brine  in  sea  ice  is  proportional  to  the 
salinity  of  the  ice  and  may  be  computed  from  the  numerical  formulas 
as  given  by  Frankenstein  and  Garner  [26]  in  parts  per  thousand: 


Vb  =  S(-52.56/T  -  2.28) 

-0.5°  >  T  >  -2.06°  C 

Vb  =  S ( -45 . 91 7/T  +  .930) 

-2.06°  >  T  >_  -8.2°  C 

Vb  =  S ( -43 . 795/T  +  1.189) 

-8.2°  _>T  >  -22.9°  C 

(2-1) 

For  temperatures  close  to  the  melting  point  an  equation  was  ob¬ 
tained  from  the  data  given  in  Assur  (1960): 

Vfa  =  SO.0028/T2  -  48. 6837/T  +  4.09201  )  -0.1°  >  T  >_  -0.5°  C  (2-2) 
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The  temperature  dependence  of  the  relative  volume  of  brine  liquid 
is  shown  in  Figure  2.2-3.  At  temperatures  above  -23°  C  sea  ice  con¬ 
sists  primarily  of  fresh  water  ice  and  brine  liquid.  For  example, 
at -10°C  l°/oo*  sea  ice  contains  993.4  grams  of  fresh  water  ice  and  6.6 
grams  of  brine  liquid.  At  -2°  C  the  brine  liquid  content  has  risen  to 
26.6  grams. 

An  empirical  equation  for  the  freezing  temperature  of  sea  water  is 
given  by  Kester  [24]: 

Tf  =  -0.137  -  .05199  S$w  -  .00007225  S^w  (2-3) 

The  freezing  point  is  important  in  the  calculation  of  the  temperature 
of  the  ice  at  the  ice-water  interface. 

2.3  Electrical  Properties 

The  electrical  properties  of  sea  ice  are  an  important  link,  along 
with  surface  roughness  and  internal  structure  of  the  medium,  in  the 
understanding  of  the  scattering  of  radar  from  it.  To  completely  de¬ 
scribe  the  radar  backscatter  properties  of  sea  ice,  at  least  a  four- 
layer  model  may  be  necessary  in  the  description  (see  Figure  2.3-1). 

The  inclusion  of  the  snow  layer  may  be  especially  important  during  the 
warmer  summer  months.  The  electrical  properties  of  sea  ice  are  directly 
linked  to  its  physical  properties.  It  is  then  evident  that  as  the 
physical  properties  change  with  time  and  age,  the  electrical  properties 
are  also  evolving. 


★ 

7„  =  parts  per  thousand. 
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A  rudimentary  four- layer 
scattering  model  of  sea  Ice 
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The  model  which  describes  the  entire  sea  ice  medium  includes  a 
snow  layer,  fresh  water  ice,  brine  liquid,  and  sea  water.  The 
dielectric  properties  at  microwave  frequencies  of  each  of  these  will 
be  reviewed. 

2.3.1  Snow 

In  most  cases  it  has  been  found  that  when  the  relatively  shallow  Arctic 
snow  layer  is  dry  it  is  transparent  to  microwave  radiation.  Prelimin¬ 
ary  results  obtained  during  the  AIDJEX  1975  experiment  indicate  that 
vertically  polarized  microwave  measurements  are  nearly  independent  of 
the  amount  of  dry  snowpack  (at  -13°  C)  encountered.  Horizontally 
polarized  microwave  measurements  were  found  to  be  strongly  effected. 

The  effect  of  liquid  water  in  the  snow  layer  has  a  strong  effect  on 
the  microwave  emission  from  the  snowpack. 

The  dielectric  constant  of  dry  snow  may  be  described  by  the  Weiner 
theory  for  mixtures  [25],  which  relates  the  dielectric  constant  to  the 
density  of  the  dry  snow.  The  relationship  between  dielectric  constant 
and  loss  tangent  of  snow  to  snow  density  is  illustrated  in  Figure  2.3-2 
by  Evans  ["26]  for  the  data  obtained  from  several  experimenters. 

Cummings  [27]  found  that  the  real  part  of  the  dielectric  constant  for 
dry  snow  was  temperature  independent  with  the  loss  tangent  sensitive 
to  temperature  as  described  in  Figure  2.3-3.  Skin  depths  calculated 
from  dielectric  constant  information,  both  experimental  and  theoretical, 
as  reported  in  the  literature  are  shown  in  Table  2.3-1. 

2.3.2  Fresh  Water  Ice 

Sea  ice  consists  of  fresh  water  ice  and  small  brine  inclusions.  The 
complex  permittivity  of  sea  ice  can  be  calculated  by  considering  the 
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effects  of  each  of  these  constituents  separately.  Much  work  has  gone 
into  the  investigation  of  the  dielectric  properties  of  pure  and  fresh 
water  ice  [26,  27,  29,  30,  31].  At  frequencies  below  300  MHz  the 
permittivity  follows  closely  the  Debye  dielectric  equation  for  polar 
molecules  and  single  relaxation  time  [26]: 

£P  =  £p  -  j£p  =  £poo  =  £poo  +  (£ps  -  V)/{1  +  jWV  (2‘4) 


where: 

Ep  =  permittivity  of  pure  ice 

Ep  =  dielectric  constant  of  pure  ice 

£p  =  loss  factor  of  pure  ice 

ep«  =  high  frequency  dielectric  constant  of  pure  ice 
£pS  =  static  dielectric  constant  of  pure  ice 
Tp  =  relaxation  time  of  pure  ice 
w  =  angular  frequency 

The  high  frequency  dielectric  constant  has  been  observed  at  microwave 
frequencies.  It  is  found  to  be  on  the  order  of  3.15  and  independent 
of  temperature.  The  loss  factor  is  mildly  temperature  sensitive,  as 
shown  in  Figure  2.3-4. 

2.3.3  Brine  Liquid 

Brine  liquid  is  a  concentrated  form  of  sea  water.  There  are  more  than 
44  elements  with  6  major  ions  (sodium,  magnesium,  potassium,  calcium, 
sulphate,  and  chloride)  in  its  composition.  Salts  increase  the 
dielectric  loss  of  water  by  adding  free  charge  carriers.  The  loss 
factor  of  sea  water  includes  an  additional  term  due  to  the  dc  conduc¬ 
tivity.  The  dielectric  properties  of  sea  water  follows  a  modified 
Debye  equation  [32]: 
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esw  =  £s=°  +  (esws  •  £s«)/(1  +  JWTsw}  *  j  asw  7  £ow  (2‘ 

where : 

e$w  =  permittivity  of  sea  water 

£sw=°  =  frequenc^  dielectric  constant  of  sea  water 
£$ws  =  static  dielectric  constant  of  sea  water 
tcw  =  relaxation  time  of  sea  water 
a...  =  dc  conductivity  of  sea  water 

5W 

eQ  =  dielectric  constant  of  free  space 

The  parameters  of  the  modified  Debye  equation  can  be  represented  as 
functions  of  temperature  only  as  illustrated  in  Figures  2.3-5  to 
2.3-9. 

2.3.4  Sea  Ice 

Measurements  of  the  dielectric  properties  of  sea  ice  and  NaCl  ice  at 
microwave  frequencies  have  been  reported  in  the  literature.  Hoekstra 
and  Cappillino  [33]  measure  the  loss  factor  of  ice  samples  obtained 
by  flash  freezing  sea  water  at  400  MHz,  9.8  GHz,  and  23  GHz.  Salin¬ 
ities  ranged  from  4  °/00  to  12  °/00  and  temperatures  ranging  from 
-60°  C  to  -30°  C  up  to  about  -5°  C.  Byrd  [34]  reported  measurements 
at  26-40  GHz  at  2  perpendicular  directions.  Temperatures  ranged  from 
-7°  C  to  -32°  C  and  salinities  ranged  from  2.85  °/oo  to  7.2  °/00. 

Bogorodsky  and  Khokhlov  [35]  made  measurements  of  the  Bering  Sea  ice 
at  10  GHz.  Salinities  ranged  from  2.5  °/00  to  7.5  °/00  and  tempera¬ 
tures  from  -13°  C  to  -10°  C. 

Hoekstra  and  Cappillino  [33],  Vant  et  al .  [36],  and  Hallikainen 
[37]  made  measurements  investigating  the  permittivity  and  loss  factor 
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6.  Relaxation  time  of  brine  and  water  as  a  function  of 
temperature. 


Conductivity  of  brine  as  a  function  of  temperature 
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8.  Dielectric  constant  of  brine  and  water  as  a  function  of 
temperature. 


of  NaCI  ice  at  microwave  frequencies  over  a  variety  of  salinities  and 
temperature  ranges. 

Hoekstra  and  Cappillino  [33]  obtained  the  following  loss  factor 
for  sea  ice  as  a  function  of  temperature  (see  Figure  2.3-10)  and  a 
function  of  frequency  (see  Figure  2.3-11).  It  is  compared  with  the 
theoretical  loss  factor  calculated  using  a  random  oriented  elongated 
brine  pocket  model  [38]: 


:si/£p  =  Ki  +^si>-/£p  =  1  -  JVbr 


*  -b/ep3 


[1  +  4esi/(esi  +  £b)] 


where: 

e  •  =  permittivity  of  sea  ice 

ep  =  permittivity  of  pure  ice 

£b-  =  permittivity  of  brine  liquid 

Vbr  -  relative  volume  of  brine  liquid  in  sea  ice 

as  a  function  of  brine  volume  in  ice  has  been  studied  by  Hoekstra  and 
Cappillino  [33].  They  measured  the  dielectric  constant  over  a  wide 
range  of  salinities  and  temperature  at  9.8  GHz.  Their  results,  as  shewn 
in  Figure  2.3-13,  draw  the  conclusion  that  dielectric  constant  is  closely 
related  to  brine  volume. 

The  theoretical  loss  factor  using  the  random  orientation  model  at 
4.7  GHz  is  shown  in  Figure  2.3-14  as  a  function  of  temperature  and 
salinity.  The  theoretical  dielectric  constant  of  sea  ice  is  shown  as 
a  function  of  temperature  and  salinity  in  Figure  2.3-15.  Figure  2.3-16 
shows  the  theoretical  attenuation  factor  of  sea  ice  at  4.7  GHz. 


Figure  2.3-12.  Theoretical  dielectric  constant  of  sea  ice  compared  to 
experimental  results  by  Vant  et  al.  (1974).  The  fre¬ 
quency  is  10  GHz. 


Figure  2.3-13.  Dielectric  constant  of  sea  ice  as  a  function  of  brine 

volume  at  9.8  GHz.  The  theoretical  results  were  obtained 
by  using  salinities  from  4  °/00  to  8  °/oc  and  temperatures 
from  -3°C  to  -16°C. 


Figure  2.3-15.  Theoretical  dielectric  constant  of  sea  ice  at  4.7  GHz  as 
a  function  of  temperature  and  salinity  (random  orientation 
model ) . 


Figure  2.3-15.  Theoretical  attenuation  factor  of  sea  ice  at  4.7  GHz  as 

a  function  of  temperature  and  salinity  (random  orientation 
model ) . 


3.0  INVESTIGATION  OF  SEA  ICE 

3.1  Introduction 

Previous  programs  to  measure  radar  backscatter  from  sea  ice  have 
involved  radar  scatterometers  flown  in  aircraft,  with  support  from 
sparse  surface-based  ice  measurements.  Numerous  imaging-radar  flights 
have  also  been  made,  but  all  of  the  systems  were  uncalibrated.  The 
aircraft  systems  used  have  been  those  that  were  available  for  use  which 
severely  restricted  attempts  to  determine  optimum  frequency  and  polar¬ 
ization.  For  these  reasons,  surface-based  scatterometer  experiments 
were  proposed  by  the  University  of  Kansas  to  provide  well -controlled 
systematic  studies  to  relate  radar  backscatter  return  to  sea  ice  and 
to  pin  down  some  of  the  uncertainties  associated  with  the  results  from 
the  various  airborne  experiments. 

3.2  May  1977  and  April  1978  Surface-Based  Missions 

The  need  to  design  an  optimum  ice  reconnaissance  radar  system  has 
been  wel 1 -documented  by  both  civilian  and  military  users.  To  do  this 
the  effects  of  polarization,  frequency,  incidence  angle,  and  resolution 
must  be  described  in  terms  of  their  contribution  in  the  ability  to  use 
radar  to  discriminate  among  types  of  sea  ice  and  to  measure  other 
properties  of  the  ice.  During  the  months  of  May  1977  and  April  1978, 
missions  were  conducted  to  contribute  to  this  description  [40,  41]. 
Active-microwave  backscatter  measurements  were  made  on  shorefast  sea 
ice  off  the  coast  of  Point  Barrow,  Alaska,  and  on  two  frozen  fresh-water 
lakes  near  the  coast.  The  objectives  of  these  measurement  programs 
were  to  study  the  available  categories  of  sea  ice,  using  in  the  study 
the  radar  parameters  of  frequency,  transmit-receive  polarization,  and 
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incidence  angle,  with  simultaneous  collection  of  ground-truth  parameters. 
Backscatter  measurements  were  conducted  using  a  four-antenna  frequency- 
modulated  continuous-wave  1-2  GHz  and  8-18  GHz  microwave  active  spec- 
trometer/scatterometer  mounted  on  a  transportable  surface-based  support 
structure  which  placed  the  antennas  at  fixed  ranges  from  the  ice. 

3.2.1  Microwave  Remote  Sensor  --  TRAMAS 
The  microwave  active  spectrometer  (MAS)  used  was  a  transportable  sur¬ 
face-based  wide-band  FM-CW  radar.  The  assembled  system  is  shown  pre¬ 
pared  for  operation  in  Figure  3.2-1.  The  L-band  (1-2  GHz)  radar  is 
located  at  the  apex  of  the  inverted  vee.  The  horn  antenna  transmits 
and  the  dish  antenna  receives.  The  Ku-X-band  (8-18  GHz)  parabolic  dish 
antennas  are  located  at  the  base  of  the  structure  and  are  aimed  at  a 
flat  plate  reflector  which  directs  the  radar  signal  to  a  focal  point 
in  the  intersection  of  the  plane  of  the  inverted  vee  and  the  ground 
plane.  Nominal  system  specifications  are  detailed  in  Table  3.2-1  and 
a  system  block  diagram  is  given  in  Figure  3.2-2. 

The  L-band  radar  was  swept  from  1.1  to  1.9  GHz  for  VV  and  VH 
polarizations.  Ku-X-band  frequencies  were  centered  at  9.0,  10.0,..., 

17.0  GHz  and  swept  +  .5  GHz.  Transmit-receive  configurations  were  HH, 

HV,  VH,  and  VV. 

Calibration  of  the  scatterometer  was  achieved  by  frequently  passing 
the  signal  through  a  delay  line  of  known  attenuation  and  by  less  fre¬ 
quently  observing  a  standard  radar  target  of  known  radar  cross-section. 

A  corner  reflector  was  used  in  the  calibration  of  the  L-band  radar  and 
a  Luneberg-lens  reflector  in  the  calibration  of  the  Ku-X-band  radar. 

Reduction  in  the  fading  effects  of  constructive  and  destructive 
interference  of  the  signal  components  of  the  scattering  elements  within 
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TABLE  3.2-1 

NOMINAL  SYSTEM  SPECIFICATIONS 


Ku-X-band 

L-band 

Type 

FM-CW 

FM-CW 

Frequency  Range 

8-18  GHz 

1.5  GHz 

Modulating  Waveform 

Triangular 

Triangular 

FM  Sweep :  Af 

1  GHz 

800  MHz 

Transmitter  Power 

14-19  dBm 

19  dBm 

Intermediate  Frequency 

50  kHz 

50  kHz 

IF  Bandwidth 

10  kHz 

10  kHz 

Antennas 

Receive  Type 

46  cm.  Reflector 

91  cm.  Reflector 

Transmit  Type 

31  cm.  Reflector 

Standard  Gain  Horn 

Feeds 

Oual  Ridge  Horn 

Log  Periodic 

Polarization 

Capabilities 

HH,  HV,  VV,  VH 

VV,  VH 

Target  Distance 

10.7  meters 

6. 1  meters 

Transmit  Beamwidth 

8.2°  at  8  GHz 

4.0°  at  17.7  GHz 

27° 

Receive  Beamwidth 

5.3°  at  8  GHz 

2.3°  at  17.8  GHz 

9.5° 

Incidence  Angle  Range 

►— * 
O 

o 

1 

o 

o 

10°  -  70° 

Cal i brat  ion 

Internal  Signal  Injection  Signal  Injection 

(delay  line)  (shorted  delay  line) 

External  Luneberg  Sphere  Square  Trihedral 

Reflector  Corner  Reflector 

-50°C  to  +50°C  -50°C  to  +50°C 


Operating  Temperature 
Range 


diagram  of  TRAMAS  system 
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the  resolution  cell  was  accomplished  by  increasing  the  number  of  inde¬ 
pendent  samples  through  frequency,  volume,  and  spatial  averaging  [42, 
43].  A  single  measurement  of  the  backscatter  from  a  particular  target 
is  almost  meaningless.  Several  samples,  preferably  a  relatively  urge 
number,  must  be  measured  and  the  results  averaged  together  to  reduce 
the  effects  of  the  fading.  The  precision  of  the  radar  scattering 
coefficient  measurement  with  a  scatterometer  is  limited  by  these  fading 
effects.  Frequency  averaging  consisted  of  sweepinq  with  bandwidth  in 
excess  of  that  required  for  resolution  [44],  Additional  independent 
looks  were  obtained  through  physically  repositioning  the  structure. 

3.2.2  Experiment  Description 

The  experiment  team  was  ground-based  out  of  the  Naval  Arctic  Research 
Laboratory  (NARL)  located  outside  of  Point  Barrow.  This  facility  pro¬ 
vided  for  messing,  billeting,  laboratory  space,  and  logistics  support. 
Backscatter  measurements  were  made  in  the  shorefast  ice  off  the  coast 
near  NARL.  Measurements  were  also  made  of  two  inland  fresh-water  lakes 
located  near  the  camp.  The  TRAMAS  system  and  support  equipment  were 
transported  to  ice  sites  by  snowmobile  and  sled  or  by  truck. 

In  1977  the  selection  of  ice  sites  was  done  with  the  aid  of  Dr. 
W.F.  Weeks,  of  U.S.  Army  Cold  Regions  Research  and  Engineering  Labora¬ 
tory  (CRREL),  by  both  ground  and  air  visual  inspection  of  the  off¬ 
shore  coastal  area  near  MARL.  Ice  categories  such  as  multiyear  and 
thick  first-year,  and  pressure  ridges,  were  located  in  the  shorefast 
ice  zone.  Inland  ice  sites  within  reach  included  fresh-water  lakes 
with  scientifically  interesting  ice  structures.  Representative  samples 
of  the  available  types  of  ice  that  were  accessible  by  snowmobile  were 


sel ected. 


59. 


Ice  conditions  in  April  1978  fell  short  of  the  ideal  conditions 
found  in  May  of  the  previous  year.  Much  of  the  shorefast  ice  zone  was 
composed  of  heavily  deformed  thick  first-year  ice.  The  site  that  was 
selected  for  the  investigation  was  a  large,  very  smooth  thick  first- 
year  ice  area  which  formed  due  to  the  refreezing  of  a  lead.  This  site 
received  much  attention  from  the  numerous  experimenters  based  at  NARL. 

3.3  Ground  Truth  Information 

The  task  of  the  ground-truth  phase  of  the  experiment  was  to  make 
measurements  at  the  time  of  the  backscatter  measurements  that  would 
describe  snow  and  ice  conditions  in  the  area  of  the  radar  footprint. 

The  following  parameters  were  used  in  the  description  of  the  snow 
pack:  depth;  density;  snow  surface  temperature;  a  general  description 
of  grain  size,  shape,  and  texture;  and  surface  roughness.  In  describing 
the  nature  of  the  sea  ice  being  studied,  the  following  parameters  were 
used:  ice  type,  thickness,  surface  roughness  (both  small  and  large 
scale),  small-scale  horizontal  inhomogeneities  (based  on  thin  sections 
from  ice  cores),  ice  surface  temperature,  ice  salinity  profile,  and 
vertical  inhomogeneities  (based  on  a  visual  core  examination).  At  two 
sites  statistical  studies  were  also  made  of  the  degree  of  preferred 
crystal  orientation  at  various  levels  in  the  ice  sheet. 


4.0  DATA  ANALYSIS 


4 . 1  Introduction 

The  normalized  measure  of  backscatter  return  is  described  by  a 
scattering  coefficient  (scattering  cross-section  per  unit  area  or  also 
called  the  differential  scattering  cross-section)  o.  This  is  a  measure 
of  the  strength  of  the  observed  signal  which  evolves  due  to  the  scatter¬ 
ing  properties  of  the  observed  region  [11].  A  discussion  of  the  scat¬ 
tering  equation  may  be  found  in  the  following  section.  The  angular 
and  frequency  responses  of  the  measured  ice  categories  are  also  presented 
and  discussed.  Comparisons  are  made  to  describe  the  ability  of  radar 
to  discriminate  among  ice  types.  The  effects  of  snow  depth,  temperature, 
and  anisotropic  properties  are  also  considered.  The  results  of  the  1-2 
GHz  measurements  are  described  in  Section  4.4  and  the  8-18  GHz  measure¬ 
ments  are  described  in  Section  4.5.  The  frequency  response  of  the  differ 
ential  scattering  cross-sections  ( c°  )  in  the  1-2  and  8-18  GHz  frequency 
regions  is  discussed  in  Section  4.6.  Application  to  ice  type/thickness 
identification  is  found  in  Section  4.7. 

4.2  Scatterometer  Equation 

The  scattering  coefficient,  c°  ,  is  used  to  describe  the  back¬ 
scatter  return  because  it  is  indeDendent  of  the  geometric  radar  oara- 
meters  (factors  which  shape  the  illuminated  area).  This  representation 
assumes  the  incoherent  collection  of  signals  from  a  very  large  number 
of  separate  scatterers  and  that  locations  of  these  scatterers  are 
sufficiently  random  in  the  resolution  cell  that  the  received  Dower  from 
each  may  be  added  to  the  power  from  all  the  others  without  consideration 

60 


4 


-1  ' 


61. 


for  phase  [45].  When  it  is  applied  in  ground-sensing  a  target  where 
the  radar  parameters  remain  essentially  unchanged  from  one  part  of  the 
resolution  cell  to  the  other,  the  radar  equation  is  given  by  the  follow¬ 
ing  [46]: 


PR  ‘ 


Pt  Gt  Gr  *2  a°  AILL 


(4tt)3 


(4 


where 


=  average  power  returned 


Pt  =  transmitted  power 

Gt  =  transmit  antenna  gain 

Gr  =  receive  antenna  gain 

A  =  wavelength 

o°  =  differential  scattering  cross-section 

A i ll  =  illuminated  area 


=  range  to  target 


Looking  at  the  average  return  power  on  a  single  scatterer  basis, 
we  have  for  a  region  which  contains  n  scattering  elements  and  which  is 
illuminated  at  one  instant  of  time  by  a  radar  the  following  equation: 
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Using  the  modified  radar  equation  (4-1),  we  may  express  the 
differential  scatterina  cross-section  a°  as: 
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The  TRAMAS  system  uses  square- law  detection  to  measure  the  return 
power.  To  calibrate  the  system,  the  returns  from  a  standard  target  of 
known  radar  cross-section  are  measured.  Immediately  before  and  after 
recording  the  return  from  the  target  of  interest,  a  coaxial  delay  line 
of  known  loss  is  switched  in  place  of  the  antenna  paths  to  provide  an 
instantaneous  measure  of  calibration.  Using  both  the  frequent  delay¬ 
line  calibrations  and  the  less-frequent  standard  radar  target  calibra¬ 
tions,  the  absolute  value  of  the  radar  cross-section  may  be  determined 
[47]. 

Numerous  return-power  measurements  of  each  ice  category  are  made. 
These  measurements  are  then  used  in  the  overall  description  of  the  ice 
category.  An  average  is  made  of  the  n  measurements  of  power  (not  dB). 
That  is,  the  n  measured  output  voltages  (v^,  v^,  v^,...,  v  )  are 
averaged  and  the  mean  a 0  is  given  (in  dB)  by: 


(dB)  =  20  log10 


1/n 


(4- 


where  v^(e)  is  the  ith  output  voltage  for  angle  of  incidence  e. 
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4.3  Surface  Truth 

The  TRAMAS  surface-based  scatterometer  is  very  attractive  because 
it  permits  easy  measurement  of  "surface  truth"  for  the  target  of  inter¬ 
est.  The  TRAMAS  system  maintains  a  fixed  range  to  target  at  all  angles 
of  incidence;  therefore,  only  the  size  of  the  resolution  cell  changes 
with  change  in  incidence  angle,  not  both  the  location  of  the  cell  and 
the  size  of  the  cell . 

The  purpose  of  the  surface-truth  task  is  to  describe  the  conditions 
of  the  target  at  the  time  of  measurement.  This  information  allows  for 
the  correlation  of  changes  in  the  physical  properties  of  the  medium  with 
changes  in  backscatter  response.  The  surface  truth  may  give  insight  on 
the  backscatter  mechanisms  and  is  very  important  if  one  compares  the 
backscatter  responses  from  similar  situations  in  the  same  experiment, 
responses  from  similar  ice  types  from  experiment  to  experiment,  or  the 
results  from  other  experiment  groups.  There  is  more  to  surface  truth 
than  the  mere  classification  of  the  ice  type,  and  it  is  hoped  that  this 
section  vividly  demonstrates  this  point  to  future  experimenters. 

4.3.1  The  May  1977  Barrow-Experiment  Surface  Truth  Measurement. 
Four  categories  of  ice  were  investigated  in  May  of  1977.  These  included 
two  thick  first-year  shorefast  sea-ice  sites,  two  fragments  of  multiyear 
sea  ice,  one  small  pressure-ridged  sea-ice  site,  one  fresh-water  inland 
lake  with  an  underlying  layer  of  water,  and  one  fresh-water  inland  lake 
frozen  to  its  mud  bottom.  The  sea  ice  sites  were  located  off  the 
Alaskan  coast  near  Point  Barrow  and  the  Naval  Arctic  Research  Laboratory 
(see  Figure  4.3-1).  The  lake  ice  sites  were  also  located  in  this 
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vicinity.  Brief  descriptions  of  these  ice  sites  may  be  found  in  Table 
4.3.1.  More  elaborate  descriptions  on  a  site-by-site  basis  follow. 

4. 3. 1.1  Site  #1:  Thick  First-Year  Shorefast  Sea  Ice.  The  ice  was 
1.65  meter  thick  and  had  approximately  10  cm  of  snow  cover.  The  surface 
was  composed  of  2-3  m  diameter  pans.  In  general,  the  topography  was 
flat  except  for  4-cm-high  vertical  steps  at  the  edges  of  the  pancakes. 

Five  sets  of  spatially-independent  backscatter  measurements  were  made. 

This  was  possible  because  the  antenna  support  structure  was  located  at 
each  of  several  positions  about  the  site.  The  air  and  snow-ice  inter¬ 
face  temperatures  for  this  site  and  the  other  sites  may  be  found  in 
Figure  4.3-2.  Snow  densities  for  this  site  and  the  other  sites  may  be 
found  in  Table  4.3-2. 

This  site  had  a  24  cm  fine-grain  slush-like  layer.  The  c-axis  orien¬ 
tation  of  these  crystals  was  apparently  random.  The  top  10  cm  had  a  dense 
bubble  distribution,  the  next  15  cm  had  a  moderately-dense  bubble 
distribution,  and  the  remainder  was  congeal ation  ice  with  a  less-dense  dis¬ 
tribution  of  bubbles  (see  Figures  4.3-3  and  4.3-4).  From  the  horizontal 
thin  section  made  at  10  cm  depth,  Figures  4.3-5  and  4.3-6,  it  is  seen 
that  the  bubble  distribution  is  random  and  that  the  bubbles  are  very 
small  in  diameter.  The  horizontal  thin  section,  when  viewed  through  a 
polarizing  filter,  demonstrates  the  random  nature  of  the  ice  crystals 
and  that  they  are  small  in  size.  The  horizontal  thin  section  made  at 
51  cm  depth  (see  Figures  4.3-7  and  4.3-8)  contains  tiny  brine  inclusions. 
These  are  systematically  distributed  and  show  a  preferred  crystal 
orientation.  When  this  section  is  viewed  through  a  polarizing  filter, 
it  is  seen  that  the  ice-crystal  platelets  at  this  level  have  grown  in  a 
preferred  direction  and  that  they  are  much  larger  than  those  in  the  upper 
parts  of  the  ice  sneet. 


Figure  4.3-5.  Horizontal  thin  section  of  Site  77-1  taken  at  depth 
10  cm. 
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Figure  4.3-6.  Horizontal  thin  section  of  Site  77-1  taken  at  10  cm 
depth  as  viewed  through  polarizing  filter. 
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Profiles  of  the  parameters  most  important  in  controlling  the  "state" 
of  the  ice  (temperature  and  salinity)  are  shown  in  Figures  4.3-9  and 

4.3- 10.  The  salinity  profile  is  a  typical  irregular  c-shaped  curve 
characteristic  of  first-year  ice  [48,49].  The  temperature  profile  is 
nearly  linear,  again  a  situation  that  is  common  in  first-year  ice  when 
there  have  not  been  recent  major  air  temperature  changes.  The  temperature 
gradient  was  .76°  C  per  m.  The  brine-volume  profile,  which  is  determined 
by  both  the  temperature  and  the  salinity  profiles  (see  equation  2-1),  is 
shown  in  Figure  4.3-9.  The  calculated  profile  is  characteristic  of  many 
other  first-year  brine-volume  profiles  that  have  been  observed  in  that 

a  slight  increase  in  brine  volume  occurs  near  the  surface,  a  minimum 
occurs  at  some  intermediate  depth  (in  this  case  15  cm),  and  a  drastic 
increase  occurs  near  the  bottom  of  the  ice  sheet.  This  latter  increase, 
of  course,  is  the  result  of  the  increase  in  salinity  near  the  interface 
with  sea  water  and  near  melting  temperatures. 

4. 3. 1.2  Site  *2:  Multiyear  Sea  Ice.  This  fragment  of  multiyear 
sea  ice  was  greater  than  three  meters  thick  with  a  snow  cover  uc  1-19 
cm.  Six  spatially-independent  backscatter  measurements  were  made  on  a 
50-cm-high  melt  hummock  of  ridge  origin  which  was  about  3  meters  in 
diameter.  The  small-scale  roughness  on  the  hummock  was  usually  4  mm. 
although  occasional  1-cm-hiqh  roughness  elements  occur.  The  upper  3  cm 
of  ice  was  granular  with  grain  size  of  1  mm  (see  Figures  4.3-12  and 

4.3- 18).  The  top  55  cm  of  the  cores  made  both  at  the  base  of  the 
hummock  (meltpond)  and  at  the  top  of  the  hummock  contained  spherical 
bubbles  up  to  4  mm  in  diameter  (see  Figures  4.3-11  to  4.3-14  and  4.3-18 
to  4.3-20).  Selow  55  cm  the  surface  was  that  of  more  typical  congealed 
ice  (see  Figures  4.3-14  to  4.3-16  and  4.3-20  to  4.3-221. 
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Figure  4.3-11.  Complete  ice  core  of  Site  77-2  (meltpond) 


Figure  4.3-12.  0-25  cm  section  of  Site  77-2  (meltpond) 


Figure  4.3-19.  20-46  cm  section  of  Site  77-2  (hummock) 


Figure  4.3-20.  41-67  cm  section  of  Site  77-2  (hummock) 


Figure  4.3-21.  64-90  cm  section  of  Site  77-2  (hummock) 


Figure  4.3-22.  89-110  cm  section  of  Site  77-2  (hummock) 


The  horizontal  thin  section  shown  In  Figure  4.3-23  shows  the  size, 
distribution,  and  shape  of  the  bubbles  In  the  upper  10  cm  of  the  Ice 
core  taken  at  the  base  of  the  hummock. 

The  horizontal  thin  section  made  at  10  cm  depth  at  the  base  of  the 
hummock  Is  shown  In  Figures  4.3-23  and  4.3-24.  The  large-diameter 
bubbles  are  found  to  be  randomly  distributed.  Bubble  shapes  vary  from 
near-perfect  spheres  to  globular  globs.  When  the  section  Is  viewed  through 
a  polarizing  filter.  It  Is  seen  that  the  Ice  crystals  at  this  depth  are 
moderate  In  size.  A  vertical  thin  section  made  at  78  cm  depth  (see 
Figures  4.3-25  and  4.3-26)  shows  a  moderately  dense  distribution  of 
bubbles  and  pockets.  These  are.  In  general,  smaller  than  those  found 
In  the  upper  layers,  with  sizes  ranging  to  the  microscopic.  This  section 
contains  many  ellipsoidal  shapes,  along  with  many  smaller  spherical 
shapes.  Again,  this  distribution  appears  to  be  random.  This  section, 
when  viewed  through  a  polarizing  filter,  shows  a  large  columnar  vertical 
crystal  growth  very  much  characteristic  of  sea  ice.  A  horizontal  thin 
section  made  at  99  cm  depth  (see  Figure  4.3-28)  contains  moderate 
crystal  sizes.  A  horizontal  thin  section  at  5  cm  depth  made  from  the 
ice  core  taken  from  the  top  of  the  hummock,  contains  less  smaller  sized 
bubbles  than  the  thin  section  made  from  the  core  taken  at  the  base  of 
the  hummock.  When  It  Is  viewed  through  a  polarizing  filter,  large  hori¬ 
zontal  crystals  are  observed  to  occur.  A  horizontal  thin  section  at  140  cm 
depth  (see  Figures  4.3-31  and  4.3-32)  contains  a  well -distributed,  moderately- 
dense  irregularly-sized  population  of  pockets  and  brine  drainage  tubes. 

Crystal  growths  are  also  larger  in  size. 


Figure  4.3-23.  Horizontal  thin  section  of  Site  77-2  (meltpond)  at 
10  cm  depth. 


Figure  4.3-24.  Horizontal  thin  section  of  Site  77-2  (meltpond)  at 
10  cm  depth  as  viewed  through  polarizing  filter. 


Figure  4.3-25.  Vertical  thin  section  of  Site  77-2  (meltpond)  at 
78  cm  depth. 


Figure  4.3-26.  Vertical  thin  section  of  Site  77-2  (meltpond)  at 
78  cm  depth  as  viewed  through  polarizing  filter. 


Figure  4.3-29.  Horizontal  thin  section  of  Site  77-2  (hummock)  at 
5  cm  depth. 
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Figure  4.3-30.  Horizontal  thin  section  of  Site  77-2  (hummock)  at 
5  cm  depth  as  viewed  through  polarizing  filter. 


Figure  4.3-31.  Horizontal  thin  section  of  Site  77-2  (hummock)  at 
140  cm  depth. 


Salinity  profiles  of  cores  taken  at  the  base  of  the  hummock,  on 
top  of  the  hummock,  and  down  the  side  of  the  hummock  are  shown  in 
Figure  4.3-33.  These  figures  show  profiles  characteristic  of  multiyear 
sea  ice  where  salinities  start  near  zero  and  increase  with  depth. 

4. 3. 1.3  Site  #3:  Thick  First-Year  Shorefast  Sea  Ice.  This  thick 
first-year  sea  ice  site  was  1.37  m  thick  with  a  uniform  2  cm  snow  cover. 
The  ice  was  very  smooth,  uniform,  and  unrafted  (see  Figures  4.3-34  to 

4.3- 39).  Six  spatially-independent  backscatter  measurements  were  made. 
The  salinity  profile,  which  is  very  characteristic  of  first-year  ice, 
is  shown  in  Figure  4.3-10. 

The  upper  2  cm  of  the  ice  sheet  (see  Figure  4.3-35)  was  a  snow- 
ice  crust  layer.  A  horizontal  thin  section  at  15  cm  depth  (see  Figures 

4.3- 38  and  4.3-39)  shows  a  light  distribution  of  moderate-sized 
channels  and  pockets  with  a  heavy,  almost  uniform,  distribution  of 
microscopic  brine  inclusions,  which  show  the  possibility  of  a  preferred 
growth  direction  of  the  moderate-sized  crystals.  A  horizontal  thin 
section  at  130  cm  depth  (see  Figures  4.3-40  and  4.3-41)  shows  a  heavy 
population  of  tiny  brine  inclusions;  the  crystals  are  large  in  size, 
and  they  exhibit  a  preferred  growth  direction.  A  vertical  thin  section 
at  122  cm  depth  (see  Figures  4.3-44  and  4.3-45)  also  shows  the  existence 
of  a  dense  population  of  brine  inclusions  and  a  pronounced  crystal 
elongation  parallel  to  the  direction  of  heat  flow,  a  strong  crystal 
orientation,  and  an  increase  in  crystal  size,  as  measured  in  the  sheet 
which  is  characteristic  of  the  columnar  zone  of  sea  ice  [48]. 

Figure  4.3-46  contains  fabric  diagrams  from  two  different  levels 
(15  and  130  cm).  Even  at  the  15  cm  level,  there  was  a  pronounced 
crystal  orientation  within  the  horizontal  plane,  and  at  the  130  cm  depth 


Salinity  -  (%o) 


Figure  4.3-34.  Complete  ice  core  of  Site  77-3 


Figure  4.3-35.  0-25  cm  section  of  Site  77-3  ice  core 


Figure  4.3-40.  Horizontal  thin  section  of  Site  77-3  at  15  cm  depth 


Figure  4.3-41.  Horizontal  thin  section  of  Site  77-3  at  15  cm  depth 
as  viewed  through  polarizing  filter. 


Figure  4.3-42.  Horizontal  thin  section  of  Site  77-3  at  130  cm  depth 


Figure  4.3-44.  Vertical  thin  section  of  Site  77-3  at  122  cm  depth 


Figure  4.3-45.  Vertical  thin  section  of  Site  77-3  at  122  cm  depth 
as  viewed  through  polarizing  filter. 
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FIGURE  4.3-46  C-axis  Orientations  Measured  Offshore  from  Barrow 
Alaska,  in  the  Chukchi  Sea  (Weeks  and  Gow,  1978). 


the  preferred  orientation  was  well  developed  (circular  standard  deviation, 
Sq  *  9.7°).  The  orientation  direction  was  almost  exactly  the  same  at  both 
levels,  with  the  c-axis  aligned  parallel  to  the  coastline  [50].  This 
orientation  was  in  excellent  agreement  with  oceanographic  observations, 
which  indicate  a  strong  northeasterly  current  parallel  to  the  Chukchi  Sea 
coast  in  the  vicinity  of  Barrow  [51,52]. 

4. 3. 1.4  Site  M:  First-Year  Pressure-Ridged  Sea  Ice.  This  site 
was  a  small  2.3-meter-high  pressure  ridge  formed  in  the  fall  when  the 
ice  was  thin.  Ice  blocks  composing  the  ridge  are  10  cm  thick.  The 
width  of  the  base  of  the  ridge  sail  is  2  m,  with  a  sail  depth  of  90  cm. 
Surrounding  undeformed  first-year  ice  is  1.3  to  1.8  m  thick.  Two  sets 
of  spatially-independent  backscatter  measurements  were  made  diagonal 

to  the  ridge.  A  cross-section  of  the  pressure  ridge  showing  the  varia¬ 
tion  in  ice  depth  and  snow  cover  in  the  vicinity  of  the  ridge  is  shown 
in  Figure  4.3-47.  The  salinity  profile  of  a  core  from  the  center  of 
the  ridge  is  shown  in  Figure  4.3-48. 

4. 3. 1.5  Site  #5:  Fragment  of  Multiyear  Ice.  This  fragment  of 
multiyear  sea  ice  was  greater  than  3  meters  thick  and  had  very  little 
snow  covering.  The  ice  surface  was  flat  on  large  scale  (in  terms  of 

a  few  meters)  with  pronounced  small  scale  melt-produced  roughness  with 
height  of  .5  to  1.0  cm.  The  upper  130  cm  of  core  had  a  layered  bubble 
structure  (see  Figures  4.3-49  to  4.3-55)  and  below  130  cm  ice  was 
normal  congealed  ice  (see  Figures  4.3-56  to  4.3-58).  A  horizontal 
thin  section  taken  at  5  cm  depth  (see  Figures  4.3-59  and  4.3-60)  shows 
a  moderately-dense  randomly-distributed  population  of  bubbles  of 
moderate  size  and  of  odd  shapes.  Crystal  diameters  ranged  in  size 


from  small  to  medium.  At  100  cm  depth  the  horizontal  thin  section 


Figure  4.3-55.  111-138  cm  section  of  Site  77-5  ice  core 


Figure  4.3-56.  132-158  cm  section  of  Site  77-5  ice  core 


Figure  4.3-59.  Horizontal  thin  section  of  Site  77-5  taken  at  5  cm 
depth. 
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Figure  4.3-60.  Horizontal  thin  section  of  Site  77-5  viewedthrough 
polarizing  filter. 


(see  Figures  4.3-61  and  4.3-62)  contains  a  lighter  population  of 
randomly-distributed  bubbles  and  brine  shafts  with  sizes  ranging  from 
very  small  (microscopic)  to  moderately  large  (in  mm's).  Ice  crystal 
sizes  were  noticeably  increased.  At  the  170  cm  depth  only  a 
moderately-dense  population  of  microscopic  brine  inclusions  could  be 
seen  in  the  horizontal  thin  section  (see  Figures  4.3-63  and  4.3-64). 

At  this  depth,  ice  crystals  were  found  to  have  grown  much  larger  than 
at  levels  closer  to  the  air-ice  interface. 

A  vertical  thin  section  was  also  made  at  this  depth  (see  Figures 
4.3-65  and  4.3-66).  In  it  the  ice  crystals  were  extremely  elongated 
and  parallel  to  the  vertical  (the  direction  of  heat  flow)  a  character¬ 
istic  of  the  columnar  zone  of  sea  ice.  Brine  drainage  tubes  were 
also  found,  which  made  themselves  apparent  through  vertical  strings 
of  small  crystals. 

The  salinity  profile  obtained  for  this  site  (see  Figure  4.3-33) 
had  the  trend  characteristic  of  multiyear  ice;  that  is,  low  salinity 
near  the  air-ice  interface,  and  an  increase  with  depth.  Five  spatially- 
independent  backscatter  measurements  were  made  at  this  site. 

4. 3. 1.6  Site  #6:  Lake  Ice  with  Ice-Water  Interface.  This  lake- 
ice  site,  Emikpuk  Lake,  had  an  ice  surface  that  was  flat  in  the  large 
scale  and  in  the  small  scale  had  an  upper  surface  that  was  very  smooth 
and  free  of  snow  cover.  This  lake  ice  was  2.06  m  in  depth  and  had  an  under¬ 
lying  water  layer.  Below  90  cm  the  ice  core  hole  filled  with  water.  Air 
and  ice  temperatures  at  time  of  the  six  spatially-independent  backscatter 
measurements  are  shown  in  Figure  4.3-2.  The  upper  layer  of  ice  was  found 
to  be  very  brittle,  which  is  characteristic  of  low-salinity  ice,  and 
thus  a  large  part  of  the  core  was  badly  disced.  The  complete  core  is 


Figure  4.3-61.  Horizontal  thin  section  of  Site  77-5  taken  at  100  cm 
depth. 


Figure  4.3-62.  Horizontal  thin  section  of  Site  77-5  viewed  through 
polarizing  filter. 


Figure  4.3-63.  Horizontal  thin  section  of  site  77-5  taken  at  170  cm 
depth. 
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Figure  4.3-64.  Horizontal  thin  section  of  Site  77-5  viewed  through 
polarizing  filter. 
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Figure  4.3-65.  Vertical  thin  section  of  Site  77-5  taken  at  170  cm 
depth. 
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Figure  4.3-66.  Vertical  thin  section  of  Site  77-5  viewed  through 
polarizing  filter. 


shown  in  Figure  4.3-67.  The  upper  9  cm  of  the  ice  sheet  (see  Figure 

4.3- 68)  had  a  heavy  population  of  moderate-sized  bubbles.  Below  9  cm 
to  about  134  cm  the  population  of  bubbles  was  very  light  and  appeared 
free  of  anything  but  the  smallest  microscopic  features.  This  purity 
was  also  apparent  in  the  thin  section  made  at  the  20  cm  depth  (see 
Figures  4.3-71  and  4.3-72).  Crystals  were  moderate  to  small  in  size. 

At  depths  between  134  and  205  cm  large  drainage  tubes  became  quite 
prominent  (see  Figure  4.3-70).  The  ice  at  these  levels  was  similar  to 
the  ice  at  the  upper  levels  except  for  the  drain  tubes.  This  is  apparent 
in  the  horizontal  thin  section  obtained  at  the  198  cm  depth  (see  Figures 

4.3- 73  and  4.3-74).  The  ice  crystals  here  are  found  to  be  very  large 
in  size. 

The  salinity  profile,  which  showed  very  low  salinity  near  the 
surface  and  increases  with  increasing  depth,  is  shown  in  Figure  4.3-75. 

4. 3. 1.7  Site  #7:  Lake  Ice  Frozen  to  Its  Mud  Bottom.  This  fresh 
water  lake  ice  site,  South  Meadow  Lake,  was  .74  meters  in  depth  and  had 
a  snow  cover  of  18  cm.  The  upper  ice  surface  was  very  smooth  and  the 
ice  sheet  extended  to  its  mud  bottom.  This  lake  was  found  to  be  very 
pure;  its  salinity  profile  is  shown  in  Figure  4.3-75.  An  attempt  was 
made  to  obtain  ice  cores,  but  because  of  the  very  low  salinity  of  the 
ice,  cores  were  badly  disced.  Three  spatial ly-independent  backscatter 
measurements  were  made  at  this  site. 

4.3.2  The  April  197C  Barrow  Experiment  Surface  Truth 

The  site  of  the  1978  surface-based  experiment,  Site  78,  was  located 


approximately  200  m  off  the  coast  of  the  Chukchi  Sea  (near  the  southwest 
end  of  the  NARL  runway).  This  site  was  selected  (also  was  limited  to) 
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due  to  the  fact  that  the  ice  along  the  Chukchi  coast  was  well  broken  up, 
no  other  ice  categories  of  interest  were  available  in  the  area,  and  there 
was  an  area  approximately  300  m  in  diameter  of  apparently  uniform,  undeformed 
fast  ice.  The  level  upper  surface  was  found  to  be  very  uniform  and  clean 
(no  visible  silt  or  sand).  The  upper  portion  of  the  ice  sheet  was  composed 
of  a  fine-grain  slush-like  lay.er  that  varied  slightly  in  thickness  from 
location  to  location.  Vertical  thin  sections  of  this  sea  ice  are  shown 
in  Figure  4.3-76.  As  can  be  seen,  the  slush  layer  here  was  4  cm  thick. 

Below  4  cm  the  crystals  immediately  elongate  parallel  to  the  vertical  (the 
direction  of  heat  flow)  and  the  elongated  brine  inclusions  characteristic 
of  sea  ice  are  evident.  The  other  three  vertical  thin  sections  show  ex¬ 
treme  crystal  elongation  which  is  characteristic  of  the  columnar  zone  of 
sea  ice.  Elongated  brine  pockets  are  particularly  evident  in  section  c 
(76  to  85  cm)  as  are  brine  drainage  tubes  which  show  themselves  through 
vertical  strings  of  small  crystals  visible  on  the  right  hand  side  of  the 
same  section  and  in  two  places  in  section  d  (153  to  167  cm). 

Horizontal  thin  sections  and  fabric  diagrams  from  the  same  block 
of  ice  are  shown  in  Figure  4.3-77.  Below  the  slush  layer  no  strong  sys¬ 
tematic  variations  in  mean  grain  size  appear  with  depth.  At  all  levels 
the  platy  substructure  of  sea  ice  that  is  revealed  by  the  systematic  dis¬ 
tribution  of  brine  inclusions  is  strongly  developed  (see  particularly  the 
thin  section  from  130  cm  depth).  Even  at  25  cm  there  is  clearly  an  iden¬ 
tifiable  preferred  crystal  orientation  in  the  horizontal  plane.  The  degree 
of  crystal  alignment  becomes  increasingly  strongly  developed  with  increasing 
depth  in  the  ice  sheet,  and  approximately  parallel  to  the  coast.  This 
is  demonstrated  by  the  systematic  decrease  with  depth  in  the  ice  of  the 


(b)  34-43  cm  (d)  153.5-162  cm 

FIGURE  4.3-76  Vertical  Thin  Sections  of  Sea  Ice  from  the  Chukchi  Sea 

Near  Barrow,  Alaska  (Site  78-17).  Sections  of  centimeter 
graph  paper  are  shown  for  scale  on  sections  (c)  and  (d). 
Photographs  made  using  crossed  polaroids. 
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standard  deviation  of  the  c-axis  orientations  about  the  mean  (measured 
in  the  horizontal  plane;  see  Figure  4.3-78).  Also,  the  mean  c-axis 
orientation  at  25  cm  is  57°  (True)  and  in  general  gradually  increases  to 
84°  near  the  bottom  of  the  ice  sheet. 

Profiles  of  the  two  parameters  most  important  in  controlling  the 
"state"  of  the  ice  (temperature  and  salinity)  are  shown  in  Figure  4.3-79. 
The  salinity  profile  is  a  typical  irregular  c-shaped  curve  which  is  char¬ 
acteristic  of  first-year  ice.  The  temperature  profile  is  nearly 
linear,  a  common  situation  in  first-year  ice  when  there  have  not  been 
recent  major  air  temperature  changes.  The  temperature  gradient  was 
10.7°  C  per  meter.  The  ice  is  not  cold  enough  for  ^01*2^0  to 
precipitate.  However,  the  solid  salt  NajSO^lC^O  should  be  present 
in  the  ice  colder  than  -8.7°  C  (in  the  ice  above  104  cm). 

The  brine  volume  profile,  shown  in  Figure  4.3-79,  was  also  deter¬ 
mined  by  the  temperature  and  salinity  profiles.  The  calculated  profile 
is  characteristic  of  many  other  first-year  brine  profiles  that  have 
been  observed  in  that  a  slight  increase  in  brine  volume  occurs  near  the 
top  of  the  ice  sheet,  a  minimum  occurs  at  some  intermediate  depth  (in 
this  case  65  cm),  and  a  drastic  increase  occurs  near  the  bottom  of 
the  ice  sheet  where  both  an  increase  in  salinity  and  near  melting 
temperatures  are  found. 

Additional  salinity  profiles  for  Site  78  are  shown  in  Figure  4.3-80. 
Photographs  of  sections  of  an  ice  core  taken  at  Site  78  are  shown  in 
Figures  4.3-81  to  4.3-87.  The  top  13  cm  of  the  ice  sheet  is  shown  in 
Figure  4.3-81.  Twelve  to  36  cm  ice  is  represented  in  Figure  4.3-82,  36 
to  85  cm  ice  is  represented  in  Figures  4.3-83  and  4.3-84,  86  to  146  cm 
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Standard  Deviation, 


Mean  C-Axis  Orientation  in  the  Horizontal  Plane  and  the  Standard 
Deviation  of  the  C-Axes  (As  Measured  in  the  Horizontal)  as  a  Function 
of  the  Vertical  Distance  Measured  from  the  Upper  Surface  of  the  Sea  Ice 


FIGURE  4.3-79  The  Salinity,  Ice  Temperature  and  Brine  Volume  Profiles  for 
the  Sea  Ice  from  Site  78-17  on  8  April  1978. 
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Figure  4.3-80.  Salinity  profiles  for  Site  78-1  and  78-2,  thick 
first-year  shorefast  sea  ice. 
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Figure  4.3-81.  0  to  13  cm  section  of  Site  78  ice  core. 


Figure  4.3-82.  12  to  25  cm  section  of  Site  78  ice  core 


Figure  4.3-84.  72  to  85  an  section  of  Site  78  ice  core 
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ice  is  represented  in  Figure  4.3-85  and  146  to  175  cm  ice  is  represented 
in  Figure  4.3-86.  The  TRAMAS  system  is  shown  studying  bare,  shorefast 
ice  in  Figure  4.3-87  and  shorefast  ice  with  a  heavy  manmade  snow  cover  in 
cigure  4.3-88. 

4.4  L-Band  Data  Analysis 

4.4.1  Results 

4. 4. 1.1  Thick  First-Year  Ice.  We  initially  believed  that  the  re¬ 
turns  from  apparently-featureless  thick  first-year  ice  should  yield 
simple  returns.  However,  as  is  shown  in  Figures  4.4-1  and  4.4-2,  this 
was  not  the  case  in  1977.  We  actually  observed  a  pronounced  clustering 
of  the  observed  returns  into  families.  This  separation  is  most  apparent 
at  large  incidence  angles.  Like-polarized  returns  at  70°  for  the  two 
groups  at  Site  77-1  and  Site  77-3  were  separated  by  12  dB  and  8  dB, 
respectively. 

The  scattering  coefficient  of  thick  first-year  ice  with  a  2  cm 
snow  cover  is  shown  in  Figure  4.4-3.  These  returns  were  found  to  be 
nicely  behaved  with  the  mean  return  falling  to  within  1-2  dB  of  a 
linear  angular  response.  When  returns  were  compared.  Site  77-1,  Site 
77-3,  and  Site  77-8  were  found  to  fall  into  three  apparent  families 
(see  Figure  4.4-4).  Returns  from  Site  78  fall  into  the  family  with 
the  lowest  return  with  Site  77-3  (3, 4, 5, 6)  and  they  follow  very  similar 
trends,  from  angles  of  20°  and  greater.  Cross-polarization  returns 
(see  Figure  4.4-5)  have  a  less  defined  grouping  into  families,  but 
may  be  said  to  group  into  two  general  categories. 
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Figure  4.4-2.  Scatterinq  coefficient  of  thick  first-year  ice.  Site  77-3, 
at  1-2  GHz,  VV. 
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Figure  4.4-3.  Scattering  coefficient  of  thick  first-year  ice,  Site  78,  at 
1-2  GHz,  VV.  Measurements  with  2-4  cm  snow  cover  were 
selected  (9  looks). 
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Figure  4.4-4.  Scattering  coefficient  of  thick  first-year  ice;  Site  77-1, 
Site  77-3,  and  Site  78;  at  1-2  GHz,  VV. 
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Figure  4.4-5.  Scattering  coefficient  of  thick  first-year  ice;  Site  77-1, 
Site  77-2,  and  Site  78;  at  1-2  GHz,  VH. 
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An  average  of  the  lower  two  families  of  the  1977  sites  was  made 
(using  <j°  but  reporting  result  as  o^g);  it  is  used  to  represent 
thick  first-year  ice  for  a  composite  1977  site,  Site  77-TFY.  A  com¬ 
parison  of  the  average  scattering  coefficient  for  Site  77-TF ?nd  Site 
78-TFY  is  shown  in  Figure  4.4-6.  The  average  angular  responses  ha.’e 
slightly  different  trends.  Scattering  from  Site  77-TFY  has  an  expon¬ 
ential  response  with  a  slope  that  changes  from  very  large  (10-45°)  to 
slight  (45-70°),  while  Site  78-TFY  has  a  rapidly-decaying  linear  response. 

4. 4. 1.2  Multiyear  Ice.  Analysis  of  multiyear  ice  was  expected  to 
be  difficult  using  data  collected  by  a  surface-based  system  because  of 
the  characteristically  non-uniform  upper  surface  of  such  ice.  The 
scattering  coefficients  of  multiyear  ice  are  shown  in  Figures  4.4-7 

and  4.4-8.  The  data  did  not  exhibit  the  expected  wide  variation  between 
independent  looks  and  was  found  to  cluster  more  tightly  than  would  be 
expected  from  surface  scatter  alone.  This  is  especially  apparent  in 
the  returns  from  Site  77-5.  A  comparison  of  multiyear  ice  sites  shows 
similar  angular  response  trends  with  absolute  levels  falling  within 
2  dB. 

4. 4. 1.3  Small  Thick  First-Year  Pressure  Ridge.  The  scattering 
coefficient  of  a  small  first-year  pressure  ridge  is  shown  in  Figure 
4.4-10.  Large  returns  were  measured  throughout  all  angles  of  incidence 
Returns  at  angles  less  than  40°  were  almost  constant,  while  returns  at 
angles  greater  than  40°  decayed  much  more  rapidly. 

4. 4. 1.4  Lake  Ice.  The  scattering  coefficients  of  lake  ice  with 
an  underlying  layer  of  water  and  of  lake  ice  that  was  frozen  to  its 
mud  bottom  are  shown  in  Figures  4.4-11  and  4.4-12,  respectively.  Both 
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Figure  4.4-6.  Comparison  of  average  scattering  coefficients  for  Site  77  - 
(1+3)  and  Site  78  at  1-2  GHz,  VV  and  VH. 
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Figure  4.4-7.  Scattering  coefficient  of  multiyear  ice,  Site  77-2,  at 


Radar  Cross  Section  -o°  (dB) 


0  10  20  30  40  50  60  70 


Angle  of  Incidence  (Degrees) 

Figure  4.4-8.  Scattering  coefficient  of  multiyear  ice,  Site  77-5,  at 
1-2  GHz ,  VV. 
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Figure  4.4-11.  Scattering  coefficient  of  lake  ice  with  underlying  layer 
of  water,  Site  77-6,  at  1-2  GHz,  VV. 
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Figure  4.4-12.  Scattering  coefficient  of  lake  ice  frozen  to  its  mud  bottom, 
Site  77-7,  at  1-2  GHz,  VV. 
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FIGURE  4.4-14  Scattering  Coefficient  of  Thick  First-Year 
and  Multiyear  Ice  at  1-2  GHz. 
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Figure  4.4-16.  Average  scattering  coefficient  of  thick  first-year,  multiyear, 
fresh  water  lake,  and  pressure  ridge  ice  at  1-2  GHz,  VH. 


sites  had  scattering  that. was  low  in  absolute  level  and  decayed  expon¬ 
entially.  These  returns,  like  the  returns  from  multiyear  ice,  did  not 
exhibit  the  expected  wide  variation  between  independent  looks;  they  were 
found  to  cluster  significantly  more  tightly  than  would  be  expected 
from  surface  scatter  alone.  This  is  especially  evident  at  the  near¬ 
vertical  angles.  The  scattering  from  the  two  lake  ice  sites  was 
almost  identical  from  10°  to  40°.  However,  after  40°  the  returns  from 
the  ice  that  was  frozen  to  its  mud  bottom  continue  to  fall  off,  while 
the  returns  from  the  ice  with  an  underlying  layer  of  water  show  a 
leveling-off  effect  (see  Figure  4.4-13). 

4.4.2  Comparison  of  Scattering  Coefficients. 

At  large  angles  of  incidence,  multiyear  ice  returns  decayed  more  rapidly 
than  those  of  thick  first-year  ice  (see  Figure  4.4-14).  The  ability  to 
discriminate  between  thick  first-year  ice,  multiyear  ice,  pressure  ridged 
ice,  and  lake  ice  is  illustrated  in  Figures  4.4-15  and  4.4-16.  Multiyear 
and  thick  first-year  ice  were  not  found  to  be  discriminable  between  10° 
and  60°.  Cross-polarization  provided  the  greatest  discriminatory  capa¬ 
bility  and  may  be  useable  at  higher  angles.  Sharply-cornered  pressure- 
ridged  ice  and  the  flat  smooth-surfaced  ice  of  the  pack  could  easily  be 
distinguished  from  each  other  at  all  but  the  highest  angle  of  incidence. 
Also,  lake  ice  had  a  significantly  lower  return  (6  dB)  than  the  sea  ice 
types  that  were  studied. 

4.4.3  The  Effect  of  Snow  Depth  on  Backs catter  Return. 


The  effect  on  snow  depth  on  the  backscatter  return  from  sea  ice  was  inves 
tigated  in  the  1978  shorefast  ice  experiments.  Natural  snow  depths  were 
found  to  vary  from  2-4  cm.  Depths  of  8  and  14  cm  were  artificially 


generated.  Experiments  were  also  performed  where  the  snow  layer  was 
removed,  and  a  bare  ice  surface  remained.  The  scattering  coefficients 
for  the  5  sites  (a  site  indicates  the  positioning  of  the  ground-based  struc¬ 
ture  at.  a  unique  location  in  the  selected  ice  area)  in  1978  were  averaged 
into  two  categories:  ice  with  a  snow  layer  and  ice  with  a  bare  surface. 

In  Figure  4.4-17  returns  at  1-2  GHz  frequencies  from  thick  first-year 
shorefast  sea  ice  with  2-4  cm  snow  covers  and  temperatures  ranging 
from  -3°  C  to  -10°  C  is  shown  to  be  insensitive  to  snow  cover.  However, 
when  examined  on  a  look-by-look  basis,  results  are  more  interesting 
and  even  confusing.  In  Figure  4.4-18  are  shown  the  returns  measured 
at  Site  78-1  for  ice  of  the  same  temperature  and  0,  2,  and  8  cm  snow 
layers.  The  bare  surface  produced  the  lowest  returns  while  the  thickest 
layer  of  snow  produced  the  greatest  return.  Large  differences  are 
found  at  near- vertical  angles.  Returns  from  ice  with  the  bare  surface 
are  about  2  dB  lower  than  returns  from  ice  with  a  2  cm  snow  cover  and 
1. 5-7.0  dB  below  the  returns  from  ice  with  an  8  cm  snow  cover.  Two 
looks  were  made  at  Site  78-2  at  the  same  temperature,  but  differing 
by  one  having  a  bare  surface  and  the  other  having  4  cm  of  snow  cover. 

Their  returns  are  nearly  identical  (see  Figure  4.4-19),  as  are  the 
returns  from  ice  with  a  bare  surface  and  ice  with  14  cm  of  snow  cover 
at  Site  78-3  (see  Figure  4.4-20),  returns  from  ice  with  a  bare  surface 
and  ice  with  2  cm  of  snow  cover  in  Site  78-4  (see  Figure  4.4-21),  and 
returns  from  ice  with  a  bare  surface  and  ice  with  4  cm  of  snow  cover 
in  Site  78-5  (see  Figure  4.4-22).  Two  things,  however,  stand  out  about 
the  responses  of  Sites  78-2,  78-4,  and  78-5.  Returns  near  nadir  are 
higher  for  ice  with  a  bare  surface,  and  ice  with  a  bare  surface  has 
a  greater  deviation  in  return  about  the  average  response. 
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Figure  4.4-17.  Average  a0  versus  angle  of  incidence  (e)  curve  for 
bare- ice  and  snow-covered  ice,  1.5  GHz. 
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Figure  4.4-18.  Scattering  coefficients  of  thick  first-year  ice 
Site  78-1,  with  a  0,  2,  and  8  cm  snow  layer  at 
1-2  GHz,  VV. 
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Figure  4.4-19.  Scattering  coefficient  of  thick  first-year  ice 
Site  78-2,  with  a  0  and  4  cm  snow  layer  at  1-2 
GHz,  VV. 
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Figure  4.4-20.  Scattering  coefficient  of  thick  first-year  ice. 
Site  78-3,  with  0  and  14  cm  snow  cover  at  1-2 
GHz,  VV. 
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Figure  4.4-21.  Scattering  coefficient  of  thick  first-year  ice, 
Site  78-4,  with  0  and  2  cm  snow  cover  at  1-2 
GHz,  VV. 


Figure  4.4-24.  To  illustrate  the  overall  effect  of  snow  depth 
on  radar  return  at  1-2  GHz,  VV,  the  average  dB 
return  for  12°  is  plotted  versus  snow  depth. 
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Figure  4.4-25.  The  effect  of  ice  temperature  on  the  backscatter 
return  from  thick  first-year  shore  fast  sea  ice 
at  1-2  GHz,  VV,  for  42,  51,  60,  and  70  degree 
angles  is  shown. 
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The  overall  effect  of  snow  depth  on  radar  return  at  1-2  GHz  is 
illustrated  in  Figure  4.4-23  where  the  return  (average  of  cfjg)  for 
42,  51,  60,  and  70  degree  angles  is  plotted  versus  snow  depth.  In 
general,  little  difference  is  detected  in  return  from  ice  with  a  bare 
surface  and  ice  with  up  to  3  cm  of  snow  cover.  Ice  with  snow  covers 
,  of  4  cm  and  8  cm  yielded  higher  returns  than  ice  with  0-3  cm  snow  covers. 

|  Ice  with  a  14  cm  snow  cover  gave  lower  return  than  ice  with  0-3  cm, 

I 

contrary  to  the  apparent  trend.  The  effect  of  snow  depth  on  radar 
return  at  near-nadir  angles  is  illustrated  in  Figure  4.4-24  where  the 

j 

!  return  at  12°  is  plotted  versus  snow  depth.  It  is  seen  that,  for  ice 

i  with  snow  covers  from  0-4  cm,  the  larqest  returns  came  from  ice  with  a 

i 

|  bare  surface  and  lowest  returns  from  ice  with  a  4  cm  snow  cover.  Returns 

t 

from  8  cm  and  14  cm  snow  covered  ice  do  not  continue  this  apparent 
trend. 

Backscatter  studies  of  snow-covered  ground  also  illustrate  a 
behavior  that  is  similar  to  that  found  for  snow-covered  sea  ice  [56]. 
j  Scattering  coefficients  at  angles  away  from  nadir  were  found  to  increase 

with  increasing  snow  cover.  However,  at  depths  approaching  2  m,  returns  did 
not  increase  significantly  with  increasing  depth. 

1  4.4.4  The  Effect  of  Temperature  on  Backscatter  Return. 

i  The  effect  of  ice  surface  temperature  on  the  backscatter  return  from 

thick  first-year  shorefast  sea  ice  at  1-2  GHz  with  varying  degrees  of 
snow  cover  (0,  2,  and  4  cms)  is  illustrated  in  Figure  4.4-25.  It  is 

I 

seen  that  backscatter  is  not  greatly  effected  by  changes  in  ice  surface 
temperature  of  -3°  C  to  -9°  C.  The  scattering  cross-section  of  snow- 
covered  qround  has  also  been  found  to  be  uneffected  by  changes  of  snow 
temperature  in  this  range  [56]. 


Radar  anisotropy  of  landfast  sea  ice  has  been  demonstrated  to  exist  by 
Kovacs  and  Morey  using  the  results  of  a  1 inearly-polarized  100  MHz  im¬ 
pulse  radar  [53].  They  observed  strong  bottom  reflections  when  the 
antenna  E-field  was  oriented  parallel  with  the  c-axis  of  the  ice  crys¬ 
tals  and  no  detectable  bottom  reflection  when  the  antenna  E-field  was 
oriented  perpendicular  to  the  c-axis. 

Results  of  the  University  of  Kansas  wide-band  scatterometry  experi¬ 
ments  and  the  study  of  crystal  alignments  in  the  fast  ice  performed  by 
Weeks  and  Gow  in  1977  (1978)  [50]  and  in  1978  (1979)  [54]  indicate  that 
microwa/e  frequencies  of  1-2  GHz  were  possibly  sensitive  to  the  aniso¬ 
tropy  of  the  preferred  crystal  orientation  in  the  fast  ice,  but  that 
these  effects  are  smaller  than  1  dB.  Since  short  wavelengths  do  not 
penetrate  to  the  depths  in  the  landfast  ice  where  strong  orientation 
effects  are  observed,  the  8.5  -  17.5  GHz  microwave  frequencies  were 
found  to  lack  any  strong  correlation  between  azimuthal  orientation 
and  strength  of  radar  return. 

Study  of  c-axis  alignments  in  1978  indicated  that  the  mean  c-axis 
orientation  was  57°  (true)  at  25  cm  depths  and  in  general  gradually 
increased  to  84°  near  the  bottom  of  the  ice  sheet.  As  a  reference, 
and  also  in  consideration  of  the  penetration  capability  of  the  1  -  2 
GHz  radar,  60°  was  selected  as  the  axis  running  parallel  to  the  mean 
c-axis.  Study  of  c-axis  alignments  in  1977  indicated  that  the  mean 
c-axis  orientation  was  35°.  Table  4.4-1  presents  a  summary  of  site 
information  that  was  considered  in  the  discussion  of  the  preferred 
c-axis  orientation  effect  on  radar  backscatter  return.  The  table  in¬ 
cludes  snow  depth,  ice  surface  temperature,  the  azimuthal  orientation 


Data  point  which  is  an  outiyer  on  regressed  line. 

Slope  resulting  when  outlying  data  point  is  not  included  in  regression. 
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of  the  antenna  E-field  with  respect  to  true  north,  the  mean  c-axis 
orientation  with  respect  to  true  north,  the  orientation  of  the  antenna 
E-field  with  respect  to  the  c-axis  of  the  landfast  ice,  the  average 
radar  scattering  cross-section  at  40°,  50°,  60°,  and  70°  incidence 
angles,  and  the  slope  of  the  linear  regression  which  was  used  to  demon¬ 
strate  the  radar  anisotropy  trend.  The  average  of  the  from 
different  incidence  angles  was  taken  to  orovide  for  a  quantitative 
description  of  the  angular  response  of  the  sea  ice.  The  40°,  50°,  60°, 
and  70°  angles  were  chosen,  since  their  return  resulted  from  the 
greatest  number  of  independent  samples  of  measurement  due  to  the 
large  illuminated  area  of  the  radar  footprint  at  these  angles  and  the 
effect  of  having  swept  in  frequency  with  excess  bandwidth. 

The  results  of  the  1977  experiment  are  shown  in  Figure  4.4-26. 

In  this  experiment  returns  were  found  to  cluster  into  two  families. 
Within  each  family  trends  indicate  that  antenna  E-field  orientations 
near  the  c-axis  of  the  landfast  sea  ice  result  in  only  a  slightly 
higher  backscatter  return  than  orthogonal  orientations.  The  slopes  of 
these  trends  were  -.0037  and  -.0025. 

The  results  of  the  1978  experiments  (see  Figure  4.4-27)  illustrate 
the  apparent  radar  response  trend  with  aximuthal  orientation  better 
than  those  for  the  1977  data.  This  may  be  due  to  the  greater  number 
of  measurements  made  in  1978.  These  returns  were  from  9  looks  with 
5  orientations.  Snow  cover  ranged  from  2  to  4  cm  and  temperatures 
ranged  from  -3°  to  -10°  C.  These  data  were  linearly  regressed,  and  a 
negative  slope  of  -.0130  dB/°  resulted. 

The  sample  correlation  coefficient  was  0.623.  To  test  the  hypo¬ 
thesis  that  there  is  zero  correlation  between  the  scatterinq  cross- 
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section  strength  and  orientation  with  c-axis,  a  significance  test  at 
the  5  percent  level  was  performed.  The  critical  value  of  this  test 
was  0.666.  This  hypothesis  was  not  rejected  at  the  5  percent  level, 
but  may  be  at  a  slightly  greater  level.  The  standard  error  of 
estimates,  sy/,  ,  which  is  due  to  the  scatter  of  th.  a0  about  the 
regression  line,  was  0.175  dB.  If  it  is  assumed  that  for  a°'s  various 
orientations  have  normal  distributions  about  the  true  regression  line 
with  common  standard  deviation  a  then  68  percent  of  the  combined 
populations  lie  within  a  +  ay^  of  the  true  regression  line  and  95 
percent  within  a  +  2  ay^  .  With  large  samples  we  could  expect  to 
approximate  this  situation. 

The  data  point  corresponding  with  Site  78-2  is  an  apparent  outlyer. 
If  it  is  considered  a  bad  point  and  is  not  included  in  the  regression, 
the  sample  correlation  coefficient  is  .926,  the  standard  error  of 
estimate  is  .0146  dB,  and  the  slope  is  -.0108  dB/®.  The  significance 
test  for  this  condition  has  a  critical  value  of  .707  for  the  5  percent 
level  and  a  value  of  .834  for  the  1  percent  level.  Therefore,  in  this 
case  the  hypothesis  that  there  is  zero  correlation  between  scattering 
cross-section  strength  and  c-axis  orientation  is  clearly  rejected. 

The  slope  of  the  regressed  data  with  the  outlyer  data  point 
included  would  indicate  that  a  1.17  dB  difference  would  exist  between 
aximuthal  look  directions  along  the  c-axis  and  perpendicular  to  the 
c-axis.  The  slope  of  the  regressed  data  with  the  outlyer  data  point 
excluded  would  indicate  that  a  .972  dB  difference  would  exist  between 
the  above  orthogonal  looks. 
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With  Respect  to  Ice  Crystals  C-Axis 


Figure  4.4-27.  Results  of  1978  experiment  demonstrating  the  effects 
of  a  preferred  ice  crystal  orientation  on  1-2  GHz 
like-polarized  (VV)  radar  backscatter  return. 
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4.4.6  Discussion  of  L-Band  Results 

Dielectric  properties  of  sea  ice  are  primarily  a  function  of  brine 
volume  in  the  ice,  which  in  turn  is  controlled  by  ice  temperature  and 
salinity.  Because  the  salinity  profiles  of  thick  first-year  ice  and 
Sites  77-1  an-i  77-3  were  similar,  ice  temperature  was  a  definite  index 
of  brine  volume.  Temperature  effects  on  radar  return  were  demonstrated 
as  shifts  in  return  trends  for  thick  first-year  ice.  Each  level  was 
directly  relatable  to  air  or  ice  temperature  at  time  of  backscatter 
measurement.  The  largest  returns  were  observed  when  air  and  ice  tempera¬ 
tures  were  the  highest.  Correspondingly  weak  returns  were  observed  when 
temperatures  were  at  their  lowest.  It  should  be  noted,  however,  that 
the  data  set  over  which  the  temperature  effect  is  described  is  small  and 
may  be  due  to  coincidence.  These  sites  also  have  different  histories  of 
ice  formation.  Site  77-1  was  formed  from  pancake  ice  (raised  edges 
occurred  above  the  normally  flat  surface),  had  an  average  of  11  cm  of 
snow  cover,  and  measurements  of  looks  3,  4,  and  5  were  taken  with  reason¬ 
ably  high  air  temperatures  (higher  air  temperatures  increase  the  free 
water  content  of  the  snow  layer)  which  may  account  for  the  wider  separa¬ 
tion  between  its  families,  whereas  Site  77-3  is  very  smooth  unrafted 
ice  with  a  2  cm  snow  cover. 

Average  backscatter  response  (a0  dB)  of  thick  first-year  ice  was 
found  to  decay  exponentially  with  incidence  angle  in  1977  and  linearly 
in  1978.  Average  returns  at  35°  to  42°  were  very  similar  in  absolute 
level . 

For  multiyear  ice,  for  the  pressure  ridge,  and  for  lake  ice,  the 
individual  looks  clustered  more  tightly  than  would  have  been  expected 
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from  fading  produced  by  surface  scatter  alone;  hence  we  conclude  that 
volume  scatter  plays  a  major  role  in  the  scattering  process,  thereby 
allowing  the  wide  sweep  width  to  give  non-coherent  averaging  of 
numerous  independent  samples  for  each  look. 

A  firm  conclusion  as  to  the  ability  to  detect  where  a  fresh  v/ater 
inland  lake  is  frozen  to  its  mud  bottom  may  not  be  possible  from  the 
backscatter  responses  obtained  from  Site  77-6  and  77-7.  Other  conditions 
were  significantly  different  that  strong  comparisons  may  not  be  valid. 
These  other  conditions  were  as  follows  for  sites  77-6  and  77-7:  ice 
temperatures  were  -2°  to  -0.5°  C  and  -0.5°  to  0°  C,  air  temperatures  were 
-2°  to  +2°  C  and  0°  to  +4°  C,  thickness  of  snow  cover  was  0  and  18  cm, 
respectively.  The  heavy  wet  snow  cover  of  Site  77-7  along  with  high 
temperatures  may  have  masked  a  more  characterise c  backscatter  return. 

4.5  KU-X-Band  Data  Analysis 

4.5.1  Results 

4. 5. 1.1  Thick  First-Year  Ice.  The  backscatter  from  what  appeared 
at  the  time  of  measurement  to  be  featureless  thick  first-year  shorefast 
sea  ice  was  found  to  be  much  as  expected  for  frequencies  in  the  8  to  18 
GHz  frequency  range.  As  the  scatter  plots  of  the  angular  response  of 
the  radar  cross-section  indicate  for  the  9,  13,  and  17  GHz  microwave 
frequencies  (VV)  (see  Figures  4.5-1  to  4.5-3),  the  overall  angular 
response  was  toward  a  common  microwave  signature.  As  discussed  in 
Section  4. 4. 1.1,  this  was  not  the  case  for  backscatter  measured  at 
L-Band  frequency.  At  this  lower  frequency  there  was  an  apparent 
breaking  of  the  microwave  signature  into  separate  families.  The  effect 
of  fading  was  also  demonstrated  in  the  fluctuation  in  the  scatter 


FIGURE  4.5-1 


Scatter  Plot  of  the  Angular  Response  of  the 
Radar  Cross-Section  of  Thick  First-Year  Ice, 
Site  77-1 ,  at  9  GHz,  HH. 
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FIGURE  4.5-2 


Scatter  Plot  of  the  Angular  Response  of  the 
Radar  Cross-Section  of  Thick  First-Year  Ice, 
Site  77-1  ,  at  13  GHz,  HH. 
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between  individual  looks  at  the  same  site.  This  was  illustrated  in 
Figures  4.5-1  to  4.5-3  by  linking  the  individual  data  points  of  a  given 
look  (orientation  of  antenna  structure).  The  scatter  plots  serve  to 
demonstrate  the  need  to  obtain  numerous  independent  samples  to  describe 
accurately  the  radar  backscatter  characteristics  of  a  target,  and  that 
a  single  independent  sample  is  of  little  value  in  the  description  of 
backscatter  properties. 

In  Site  77-1  an  interesting  feature  does  generate  some  curiosity. 

At  9  GHz  and  17  GHz  there  is  a  noticeable  increase  (2.5  and  8  dB, 
respectively)  in  the  scattering  cross-section  at  50°  and  40°,  respec¬ 
tively.  This  increase  may  be  explained  by  the  action  of  the  edges  of 
the  pancake  ice  that  compose  the  top  few  cm  of  the  ice  sheet;  those  may 
act  as  corner  reflectors,  due  to  the  discontinuous  surface  and  the 
short  wavelengths. 

Comparisons  of  the  angular  response  of  the  average  radar  scattering 
cross-section  (see  equation  4-4  for  definition  of  average  o°  response) 
for  HH  and  VV  like-polarizations  and  an  average  of  the  HV  and  VH  cross¬ 
polarizations  are  shown  in  Figures  4.5-4  to  4.5-6.  These  average  radar 
cross-sections  were  also  linearly  regressed  against  angle  of  incidence 
(outlyers  were  not  included  in  the  regression)  and  the  regression  line 
was  included  in  the  figures.  Results  indicate  that  VV  and  HH  scattering 
coefficients  were  very  much  the  same  in  absolute  level  and  that  the 
linear  fit  applied  well  over  most  of  the  range  of  angle  of  incidence. 
However,  at  17  GHz,  scatter  for  the  two  like-polarizations  (HH  and  VV) 
was  separated  by  4  dB,  although  the  slopes  of  the  straight  line  fits 
were  nearly  identical.  It  should  be  noted  that  small  slope  differences 
have  very  little,  if  any,  significance.  For  example,  in  this  case,  the 


FIGURE  4.5-5  Scattering  Coefficient  of  Thick  Cirst-Year 
Ice ,  Site  77-1  ,  at  13  GHz . 
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slopes  and  corresponding  intervals  of  95%  confidence  for  VV  and  HH 
polarizations  were  0.131  +  .041  and  0.196  +  .086,  respectively.  The 
scattering  coefficient  for  cross-polarization  was  found  to  be  very 
much  lower  in  absolute  level  than  those  for  like-polarizations.  This 
separation  of  level  was  16,  14  and  7  dB  at  40°  for  9,  13,  and  17  GHz, 
respectively.  This  may  indicate  that,  as  frequency  increased,  the  de¬ 
polarizing  mechanism  of  the  medium  becomes  more  significant. 

Like-polarization  (HH  and  VV)  and  averaged  cross-polarization 
returns  for  Site  77-3  are  shown  in  Figures  4.5-7  to  4.5-9.  Interval 
bars  are  also  included  to  show  the  maximum  deviation  about  the  average 
VV  return.  Deviations  from  the  average  for  HH  and  cross-polarizations 
were  found  to  behave  very  much  the  same  as  the  deviations  found  in  the 
VV  case.  The  angular  response  was  again  found  to  be  well -fitted  with 
a  straight  line  between  the  angles  of  20°  and  70°. 

Scatter  for  the  two  like-polarizations  followed  very  much  the  same 
angular  response,  but  exhibited  differences  in  absolute  levels  of  2  to 
3  dB  between  VV  and  HH.  Scatter  for  cross-polarization  was  found  to  be 
well -separated  from  the  1  ike-polarization  scatter  -  by  13,  13.5,  and 
10  dB  at  40°  for  9,  13,  and  17  GHz,  respectively.  The  cross-polariza¬ 
tion  angular  response  also  showed  a  slope  that  was  consistently  less 
than  the  slope  of  the  like-polarization  response. 

The  radar  scattering  cross-sections  of  shorefast  thick  first-year 
sea  ice  with  natural  snow  cover  measured  at  the  one  major  site  in  1978, 
Site  78,  are  shown  in  Figures  4.5-10  to  4.5-12.  The  angular  response 
of  the  scatter  was  found  to  be  reasonably  well  fit  with  a  straight  line. 
This  fit  was  usually  found  to  be  good  from  20°  to  70°,  but  in  many 
instances  it  was  a  very  good  fit  in  to  12°  (9  GHz,  13  GHz-VV,  9  GHz- 
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FIGURE  4.5-7  Scattering  Coefficient  of  Thick  First' 
Year  Ice,  Site  77-3,  at  9  GHz. 
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Radar  Cross  Section  - 


FIGURE  4.5-9  Scattering  Coefficient  of  Thick  First-Year  Ice, 
Site  77-3,  at  17  GHz. 
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FIGURE  4.5-10  Scattering  Coefficient  of  Thick  First- Year  Ice 
with  Natural  Snow  Cover,  Site  7,  at  9  GHz. 
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FIGURE  4.5-11  Scattering  Coefficient  of  Thick  First-Year  Ice 
with  Natural  Snow  Cover,  Site  78,  at  13  GHz. 
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FIGURE  4.5-12  Scattering  Coefficient  of  Thick  First-Year  Ice 
with  Natural  Snow  Cover,  Site  78,  at  17  GHz. 
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Cross,  13  GHz-HH,  and  17  GHz-Cross).  The  like-polarization  angular 
response  of  VV  and  HH  had  nearly  the  same  slopes,  but  differed  in 
absolute  level  (0.5,  1.5,  and  2.0  dB  for  9,  13  and  17  GHz,  respectively). 
Again,  cross-polarization  scatter  exhibited  slopes  that  were  slightly 
less  steep  than  the  slopes  of  like-polarization.  Separations  between 
like  (HH)  and  cross-polarization  (HV)  were  found  to  be  16,  12,  and  11 
dB  at  40°  for  9,  13,  and  17  GHz,  respectively. 

Comparisons  of  Site  77-1,  Site  77-3,  and  Site  78  at  9,  13,  and  17 
GHz  (HH)  are  shown  in  Figures  4.5-13  to  4.5-15.  The  overall  comparison 
of  the  three  angular  responses  for  the  three  thick  first-year  sea  ice 
sites  of  slightly  different  origins  shows  similar  absolute  levels  and 
similar  slopes.  Site  78  was  found  to  have  an  absolute  level  that  was 
consistently  lower  than  the  levels  found  in  the  other  two  sites.  It 
had  an  angular  response  that  was  fairly  well-behaved  (with  respect  to 
the  regressed  straight  line).  The  angular  response  of  Site  77-3  was 
also  well-behaved,  but  the  angular  response  of  Site  77-1  varied  greatly 
about  its  straight  line  fit. 

4. 5. 1.2  Small  Thick  First-Year  Pressure  Ridge.  The  angular  re¬ 
sponse  of  the  radar  cross-section  (HH,  VV,  and  Cross)  at  13  GHz  is 
shown  in  Figure  4.5-16.  Responses  for  HH  and  VV  like-polarizations  were 
very  much  the  same  with  nearly  identical  slopes,  although  separated  in 
absolute  level  by  3  dB.  Returns  at  70°  were  much  lower  than  the  returns 
at  all  other  angles  and  fell  significantly  away  from  the  trend  line. 

This  response  may  possibly  be  peculiar  to  the  geometry  of  the  measure¬ 
ment  scheme  and  the  pressure  ridge  itself. 

4. 5. 1.3  Multiyear  Ice.  The  average  angular  response  (with  devia¬ 
tions)  of  the  radar  backscatter  cross-section  of  multiyear  ice.  Site 
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FIGURE  4.5-13  A  Comparison  of  the  Scattering  Coefficients  of 
Thick  First-Year  Sea  Ice  for  Sites  77-1,  77-3, 
and  78  at  9  GHz,  HH. 
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FIGURE  4.5-14  A  Comparison  of  the  Scattering  Coefficients 
of  Thick  First-Year  Sea  Ice  for  Sites  77-1, 
77-3,  and  78  at  13  GHz,  HH. 
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FIGURE  4.5-15  A  comparison  of  the  scattering  coefficients  of  thick 
first-year  sea  ice  for  Sites  77-1,  77-3,  78  at  17  GHz, 
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FIGURE  4.5-16 


Scattering  coefficient  of  a  small  thick  first-year 
pressure  ridge,  Site  77-4,  at  13  GHz. 
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77-2  and  Site  77-5,  is  shown  for  9,  13,  and  17  GHz  (HH)  in  Figures 
4.5-17  to  4.5-19. 

The  scattering  coefficients  from  the  two  fragments  of  multiyear 
ice  overlap  reasonably  well.  This  is  interesting  because  these  two 
sites  were  examples  of  two  different  kinds  of  surface  roughness  that 
may  occur  in  multiyear  ice.  Site  77-2  was  of  a  large  hummock  formation 
(a  large  amplitude  and  prominent  rolling  feature)  while  Site  77-5  would 
be  described  as  having  rolling  features  of  small  amplitude.  The  higher 
the  frequencies  the  greater  was  the  spread  between  the  average  micro- 
wave  signatures  (1.5,  2,  and  4  dB  at  40°  and  9,  13,  and  17  GHz,  respec¬ 
tively). 

4. 5. 1.4  Lake  Ice.  Average  scatter  for  lake  ice  with  an  underlying 
layer  of  water.  Site  77-6,  is  shown  in  Figures  4.5-20  to  4.5-22.  Re¬ 
sults  were  fit  with  a  linear  regression  of  all  or  selected  data  points. 
Both  like-polarizations  exhibited  almost  identical  slopes,  but  they 
showed  separations  in  absolute  level  of  2.5,  1.5,  and  3.2  dB  at  9,  13, 
and  17  GHz,  respectively.  Again,  the  cross-polarization  response  showed 
a  smaller  slope  than  the  like-polarization  responses,  and  they  were 
separated  in  absolute  level  at  40°  by  8,  11,  and  6  dB  for  9,  13,  and 
17  GHz,  respectively.  Interestingly,  there  was  a  significant  drop  in 
return  strength  for  all  polarizations  at  70°  (1,  3,  and  2  dB  at  9,  13, 
and  17  GHz,  respectively). 

Scatter  for  lake  ice  which  was  frozen  to  its  mud  bottom.  Site  77-7, 
is  shown  in  Figures  4.5-23  to  4.5-25.  The  angular  responses  were  linear 
in  shape,  this  being  especially  true  at  angles  of  incidence  from  35°  to 
70°.  HH  and  VV  responses  were  nearly  identical,  with  separation  in 
absolute  level  being  less  than  one  dB  throughout  the  frequency  range. 
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FIGURE  4.5-18  Scattering  coefficient  of  multiyear  ice,  Site  77-2 
and  Site  77-5,  at  13  GHz,  HH. 
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FIGURE  4.5-20  Scattering  coefficient  of  lake  ice  with  underlying 
layer  of  water.  Site  77-6,  at  9  GHz. 
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FIGURE  4.5-21  Scattering  coefficient  of  lake  ice  with  underlying 
layer  of  water,  Site  77-6,  at  13  GHz. 
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FIGURE  4.5-25  Scattering  coefficient  of  lake  ice  frozen  to  its 
mud  bottom.  Site  77-7,  at  17  GHz. 
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The  cross-polarization  response  had  a  slope  that  was  smaller  than  those 
for  the  like-polarization  cases.  Cross-polarization  scatter  was  also 
well-separated  from  like-polarization  scatter  (14,  14,  and  10  dB  for 
9,  13,  and  17  GHz,  respectively). 

Comparison  of  VV  scatter  for  Site  77-6  and  Site  77-7  (see  Figures 
4.5-26  to  4.5-28)  shows  a  wide  separation  between  coefficients  for  these 
ice  categories  (5,  6,  and  7  dB  for  9,  13,  and  17  GHz,  respectively).  In 
all  cases,  the  return  from  the  lake  which  was  frozen  to  its  mud  bottom, 
which  also  had  a  heavy  snow  cover  and  higher  air  and  ice  temperatures, 
was  found  to  be  the  higher.  Response  trends  for  9  and  13  GHz  show  that 
both  sites  had  responses  that  were  nearly  identical  in  slope  but  were 
separated  in  absolute  level.  Cross-polarization  scatter  from  the  two 
sites  were  not  separated  as  dramatically  as  for  the  like-polarization 
case  (1.3,  3.5,  and  4  dB  at  40°  for  9,  13,  and  17  GHz,  respectively) . 

4.5.2  Comparison  of  Scattering  Coefficients 

Mean  scattering  coefficient  trends  of  thick  first-year  and  multi¬ 
year  ice  for  frequencies  in  the  Ku-X  band  microwave  region  are  compared 
in  Figures  4.5-29  to  4.5-31.  The  radar  parameters  are  9,  13,  and  17  GHz 
frequencies  and  HH,  HV,  and  VV  transmit-receive  polarizations.  The 
angular  response  at  9  GHz  is  shown  in  Figure  4.5-29.  The  minimum  spread 
(HH)  which  separated  returns  from  these  ice  types  was  6  dB  allowing 
them  to  be  discriminated  at  all  polarizations  and  angles.  Cross-polari¬ 
zation  provided  the  best  discrimination,  affording  12  dB  of  separation, 
which  was  independent  of  incidence  angle  from  20°  to  70°.  VV  polariza¬ 
tion  was  slightly  better  than  HH,  afforded  8  dB  of  separation,  and  was 
best  from  20°  to  45°.  The  13  GHz  response  was  similar  (Figure  4.5-30). 
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Figure  4.5-26.  Scattering  coefficierrt  of  lake  ice,  Site  77 
Site  77-7,  at  9  GHz,  VV  and  Cross. 
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Figure  4.5-28.  Scattering  coefficient  of  lake  ice,  Site  77-6 
and  Site  77-7,  at  17  GHz,  VV  and  Cross. 
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Its  angular  response  was  nearly  identical  to  the  response  at  9  GHz 
except  for  a  higher  absolute  level.  The  17  GHz  response  was  also 
similar,  but  with  less  separation  in  the  returns  from  the  two  ice  types 
(Figure  4.5-31 ) . 

Comparisons  of  the  average  scattering  coefficients  of  multiyear 
and  thick  first-year  sea  ice  (discussed  above)  with  pressure-ridged 
sea  ice  and  fresh  water  lake  ice  are  shown  vs.  angle  of  incidence  in 
Figures  4.5-32  to  4.5-37  for  9  and  17  GHz  (HH)  and  13  GHz  (HH,  VV,  and 
HV).  In  all  cases,  multiyear  ice  had  the  highest  backscatter  return 
and  would  be  easily  discriminable  from  the  other  ice  cypes.  The  pres- 
sure-ridged  ice  had  the  next  highest  backscatter  return.  It  would  be 
hard  to  discriminate  pressure-ridge  ice  from  thick  first-year  ice  at 
near  grazing  angles  on  the  grounds  of  backscatter  return  only,  but 
pressure  ridges  are  linear  features  and  thus  are  readily  identifiable  in 
radar  imagery  of  floe  features.  Lake  ice  was  found  to  have  the  lowest 
backscatter  response  of  these  ice  categories,  with  returns  well  below 
those  of  the  sea-ice  categories.  HH  and  VV  polarizations  have  similar 
responses  for  all  types  of  ice.  This  is  demonstrated  when  one  compares 
Figures  4.5-33  and  4.5-34  which  show  scatter  for  13  GHz  HH  and  VV. 
Cross-polarization  returns  are  very  low  compared  with  VV  and  HH  returns 
(Figure  4.5-35). 

From  backscatter  alone,  one  may  find  that  if  proper  angles  were 
not  chosen,  difficulty  may  arise  in  the  ability  to  discriminate  between 
pressure-ridged  and  multiyear  ice,  and  also  between  lake  ice  and  thick 
first-year  sea  ice.  However,  in  an  actual  situation,  these  problems 
would  not  impair  the  ability  to  discriminate  due  to  the  available  infor¬ 
mation  gained  through  study  of  the  tones  and  shapes  of  the  features 
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FIGURE  4.5-33  Scattering  Coefficient  of  Thick  First-Year,  Multiyear, 
Fresh  Water  Lake,  and  Pressure  Ridge  Ice  at  13  GHz. 
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under  study. 


4.5.3  The  Effect  of  Snow  Depth  on  Backscatter  Return 

Snow  was  found  to  affect  significantly  the  absolute  level  of  the 
radar  cross-section  of  thick  first-year  sea  ice.  This  effect  ranged 
from  1  to  4  dB  depending  upon  the  frequency  and  angle  of  incidence. 

Snow  was  shown  to  be  a  contributing  source  of  backscatter  from  the 
sea  ice  medium.  The  effect  of  the  snow  surface  was  to  enhance  the 
backscatter  by  modifying  the  overall  absolute  level  of  the  scatter,  with 
the  shape  of  the  angular  response  remaining  very  much  the  same.  The 
contribution  of  the  snow  layer  to  the  backscatter  response  increased 

with  frequency.  These  effects  are  illustrated  in  Figures  4.5-37  to 
4.5-39.  Similar  results  have  also  been  found  to  be  true  for  scatter 
from  snow-covered  ground  [56]. 

A  specific  example  of  the  scatter  from  ice  with  a  bare  surface  is 
shown  in  Figure  4.5-40.  This  figure  shows  the  angular  response  of 
Site  78-2  for  two  looks  at  13  GHz.  One  look  had  a  4  cm  snow  cover 
(natural)  and  the  other  had  a  bare  surface.  The  ice  with  a  bare  surface 
exhibited  an  8  dB  lower  absolute  level  than  the  ice  with  a  snow  surface. 
Two  interesting  effects  were  also  noted  in  the  response  from  a  bare- 
surface.  There  was  a  significant  increase  in  absolute  level  near  nadir 
(12°)  and  a  significant  decrease  in  return  at  the  largest  angle  (70°). 
These  are  characteristics  of  scatter  from  a  smooth  surface  and  were 
found  to  be  exhibited  in  many  data  sets.  The  overall  response  of  the 
radar  scattering  coefficient  as  it  varied  with  snow  depth  is  shown 
in  Figure  4.5-41.  In  this  figure,  the  response  of  the  radar  cross- 
section  fro  42,  51,  60  and  70°  has  been  averaged  (in  dB)  and  this 
quantitative  measure  was  used  to  serve  as  a  guide  to  indicate  the  degree 


to  which  the  response  of  backscatter  is  effected  by  snow  depth.  Results 
show  that  ice  with  a  bare  surface  exhibited  lower  scatter  than  ice  with 
a  snow  cover  (at  least  up  to  14  cm),  and  that  ice  with  snow  depths  up 
to  8  cm  showed  greater  absolute  levels  of  scatter  than  ice  with  a 
thinner  snow  cover.  The  14  cm  snow  cover  was  generated  by  shoveling 
snow  from  outside  of  the  radar  footprint  onto  the  snow  in  the  radar 
footprint,  leveling  the  piled  snow,  and  slightly  packing  the  snow  until 
the  14  cm  depth  was  obtained;  its  backscatter  was  only  slightly  higher 
than  the  scatter  measured  for  bare  ice.  Thus,  the  packing  apparently 
reduced  the  snow  backscatter;  detailed  study  of  differences  in  the  snow 
properties  was  not  made. 

4.5.4  Temperature  Effects  on  Backscatter  Return 

Measurements  were  obtained  in  1978  that  allowed  some  study  into 
the  effects  of  temperature  on  backscatter  return.  Since  the  1978  experi¬ 
ment  site  was  situated  on  one  field  of  ice  with  the  same  history  of  for¬ 
mation,  variations  in  measurements  were  correlated  with  the  accompanying 
variations  in  temperature.  This  task  was  attempted  in  1977,  but  proved 
difficult  due  to  the  limited  number  of  looks  at  a  given  site. 

Results  indicate  that,  in  general,  lower  radar  scattering  cross- 
sections  were  measured  at  lower  temperatures.  This  was  found  to  be  true 
for  ice  with  a  bare  surface  (Figure  4.5-42),  ice  with  2  cm  of  snow  cover 
(Figure  4.5-43),  and  ice  with  a  4  cm  snow  cover  (Figure  4.5-44).  Results 
are  not  conclusive  due  to  the  limited  amount  of  data  and  a  significant 
amount  of  scatter  in  this  data,  but  a  similar  response  has  been  demonstra¬ 
ted  by  snow-covered  ground  [56].  The  indication  is  that  temperature 
effects  may  be  in  the  1  to  8  dB  range  for  changes  in  temperature  of 
10°  C  depending  upon  frequency  and  snow  cover. 
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Effect  of  temperature  on  the  scattering  coefficient 
of  thick  first-year  ice  with  a  4cm  snow  cover  at  9 
and  13  GHz,  an  angle  of  incidence  of  55°,  and  HH 
polarization. 


The  frequency  responses  of  multiyear  and  thick  first-year  sea  ice 
for  like-(HH)  and  cross-(HV)  polarization  scatter  at  70°,  51°  and  22° 
angles  of  incidence  are  shown  in  Figures  4.6-1  to  4.6-4.  A  comparison 
of  the  frequency  response  of  the  scatter  at  60°  for  the  thick  first-year 
at  site  77-3  and  Site  73  is  shown  in  Figure  4.6-5.  Trends  indicate 
that  the  radar  cross-section  at  large  angles  of  incidence  increases 
linearly  with  increasing  frequency  in  the  8.5  to  17.5  GHz  region.  The 
response  at  1.5  GHz  raises  some  interesting  questions.  If  the  backscat- 
ter  mechanisms  that  apply  at  8  to  18  GHz  frequencies  were  applicable  to 
the  1  to  8  GHz  region,  one  would  only  need  to  extend  the  approximately 
linearly  increasing  frequency  response  used  to  describe  the  scatter  in 
the  8  to  18  GHz  region  into  the  1  to  8  GHz  region.  The  response  at  1.5 
GHz  shows  that  different  scattering  mechanisms  must  apply  in  the  case 
of  thick  first-year  ice.  Thick  first-year  ice  backscatter  was  higher 
at  1.5  GHz  by  greater  than  4  dB  than  the  extrapolated  value.  Multiyear 
ice  would  be  detectably  higher  than  thick  first-year  ice  by  extrapola¬ 
tion,  whereas  it  was  found  that  the  measured  1.5  GHz  return  from  thick 
first-year  ice  (Site  77-3)  was  higher  in  absolute  level  than  that  from 
multiyear  ice  (Site  77-5)  by  at  least  2  dB.  The  closest  that  multiyear 
scatter  at  1.5  GHz  comes  to  the  extended  8-18  GHz  response  trend  is 
4  dB  (measured  too  low).  Since  VV  scatter  has  been  found  to  be  consis¬ 
tently  higher  in  absolute  level  than  HH  scatter,  this  level  of  separa¬ 
tion  between  the  8  to  13  GHz  (HH)  trend  and  the  1.5  GHz  measurement  may 
be  even  more  significant.  However,  an  interesting  phenomenon  occurred 
in  the  9  to  13  GHz  regions.  Scatter  at  these  frequencies  has  a  fre¬ 
quency  response  that  is  very  linear  in  many  cases,  but  with  a  slope  that 
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FIGURE  4.6-2  Scattering  Coefficient  Frequency  Response  of  Thick  First-Year  and  Multiyear 
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Figure  4.6-5.  Scattering  coefficient  frequency  response  of  thick  first-year  sea  ice,  Site 
77-3  and  Site  78. 


is  different  than  the  slope  of  the  regression  from  8-18  GHz.  It  is 
seen  that  scatter  from  1-13  GHz  may  be  best  described  using  a  straight 
1 ine  with  this  slope. 

Cross-polarized  scatter  from  multiyear  ice  was  greater  in  absolute 
level  than  like-polarized  scatter  from  thick  first-year  ice  in  the  8  to 
18  GHz  region.  The  frequency  response  at  the  different  angles  of  in¬ 
cidence  was  found  to  be  much  the  same,  which  was  an  expected  result 
since  their  angular  responses  were  linear  in  nature. 

In  comparing  the  frequency  responses  (like-polarization)  of  the 
thick  first-year  ice  of  Site  77-3  and  Site  78  (see  Figure  4.6-5),  it 
was  found  that  they  correlated  well  both  in  frequency  variation  and 
absolute  level.  Cross-polarized  scatter  for  these  two  sites  in  two 
different  years,  however,  differed  greatly.  This  was  especially  true 
at  the  higher  frequencies. 

4.7  Application  to  Ice  Type/Thickness  Identification 

A  user  of  an  operational  ice  reconnaissance  system  needs  to  be 
able  to  use  its  output  in  describing  the  ice  distribution  of  his  area 
of  interest.  The  categories  of  ice  located  in  the  region  should  be 
identifiable  and  relatable  to  ice  thickness  ranges.  The  fraction  of 
ice  in  each  category  should  also  result  from  the  interpretation  process. 
The  radar  backscatter  data  obtained  in  the  1977  and  1973  experiments  may 
be  used  to  provide  insight  into  this  capability  for  a  wide  range  of 
frequencies,  incidence  angles,  and  polarizations.  Multiyear,  thick 
first-year,  and  pressure-ridged  sea  ice  and  fresh  water  lake  ice,  appear 
in  this  1 imited  data  set. 

At  1.5  GHz  frequency,  the  radar  did  not  have  the  ability  to  dis- 
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criminate  between  thick  first-year  and  multiyear  sea  ice  for  angles 
between  12°  and  60°;  however,  discrimination  may  be  possible  at  angles 
larger  than  60°.  With  cross  polarization  these  ice  types  could  be  more 
readily  distinguished  than  with  like  polarization  (see  Figure  4.7-1). 

The  difference  in  the  radar  cross-section  of  multiyear  and  thick  first- 
year  sea  ice  at  an  angle  of  incidence  of  70°  was  found  to  be  3  dB  for 
like-polarization  (VV)  and  5  dB  for  cross-polarization  (VH).  Returns 
at  this  frequency  allowed  easy  distinction  between  the  flat  areas  of 
the  pack  ice  and  the  more  prominent  features,  such  as  a  pressure-ridge. 
Lake  ice  returns  were  much  lower  at  least  than  the  returns  from  any  of 
the  categories  of  sea  ice  that  were  examined. 

Ku-X-Band  results  indicate  that  frequencies  between  8  and  18  GHz 
allow  discrimination  between  thick  first-year,  multiyear,  and  pressure- 
ridged  sea  ice  and  fresh  water  lake  ice.  Both  like-  and  cross-  polari¬ 
zations  were  found  to  be  useable. 

One  of  the  most  important  needs  is  the  ability  to  discriminate 
between  thick  first-year  ice  and  old  ice  (multiyear).  Old  ice  is  much 
thicker  than  first-year  ice  and  therefore  concerns  most  users  of  the 
Arctic  region.  Figures  4.7-1  and  4.7-2  graphically  display  the 
difference  in  backscatter  return  at  frequencies  of  1.5,  9,  13,  and  17 
GHz  between  these  two  ice  types.  Both  like-polarization  (VV)  and 
cross-polarization  (VH)  differences  are  shown  in  separate  graphs.  Re¬ 
sults  indicate  that  both  like-  and  cross-  polarizations  were  useable. 

The  9  GHz  frequency,  for  both  like-  and  cross-  polarizations,  was  found 
to  discriminate  best.  Multiyear  ice  returns  were  about  S  dB  higher 
for  like-polarization  (best  angles  were  20°  to  43°)  and  12  dB  higher  for 
cross-polarization  (best  angles  were  35’  to  70°)  than  the  returns  from 
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thick  first-year  ice.  Returns  at  13  GHz  and  17  GHz  were  found  to  dis¬ 
criminate  about  equally  well.  Both  were  about  2  dB  worse  than  those  at 
9  GHz.  For  cross-polarization  returns  the  difference  between  thick 
first-year  and  multiyear  ice  was  much  more  sensitive  to  frequency,  with 
difference  in  returns  for  17  GHz  being  much  less  than  those  for  9  GHz 
(6  dB). 

Another  task  of  an  operational  system  is  to  provide  information 
about  the  extent  of  pressure  ridges.  Also  desired  is  the  ability  to 
estimate  the  size  of  the  pressure  ridge  from  the  width  of  the  bright 
return  and/or  the  extent  of  the  radar  shadow  associated  with  this 
linear  feature.  The  pressure  ridge  that  was  examined  in  the  1977 
experiment  was  small,  and  the  response,  especially  at  the  large  inci¬ 
dence  angles,  may  have  been  affected  by  the  orientation  of  the  ridge 
with  respect  to  the  radar  (the  radar  looked  parallel  to  this  linear 
feature).  Pressure-ridges  typically  provide  bright  returns  on  radar 
imagery  at  large  incidence  angles  and  at  these  frequencies  the  back- 
scatter  from  the  small  ridge  fell  between  the  returns  from  multiyear 
and  thick  first-year  ice.  Indications  are  that  20°  to  40°  incidence 
angles  would  be  useful  in  locating  ridges  from  backscatter  alone. 
However,  since  ridges  are  linear  features  with  returns  greater  than 
those  of  thick  first-year  ice,  they  should  be  easily  identified. 
Separation  of  ridges  from  multiyear  ice  is  done  readily  since  multiyear 
ice  has  the  characteristics  of  a  floe  with  angular  edges  and  has  a 
bright  radar  return. 

Little  difficulty  would  exist  in  discriminating  sea  ice  from  fresh 
water  ice  (our  case  was  a  fresh  water  lake).  The  returns  for  the  lake 
that  had  very  little,  if  any,  snow  cover,  and  had  an  underlying  layer  of 


water,  were  found  to  be  at  least  6  dB  lower  than  the  returns  from  thick 
first-year  ice  (which  has  the  lowest  return  of  any  of  the  sea  ice 
categories  that  were  studied).  The  lake  that  was  frozen  to  its  mud 
bottom  and  had  a  heavy  wet  snow  cover  was  found  to  have  returns  that 
were  on  the  order  of  8  dB  higher  than  the  previously  described  lake. 

The  uncharacteristically  high  returns  appear  to  be  typical  of  ice  with 
its  snow  surface  in  a  melt  condition.  For  example,  it  has  been  found 
for  sea  ice  that  returns  may  increase  to  a  point  in  which  there  is 
little,  if  any,  ability  to  discriminate  among  types. 
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5.0  MODELING  OF  BACKSCATTER  RETURN 

Theory,  as  long  as  it  can  be  confirmed  by  experiment,  provides 
a  useful  base  from  which  one  may  judge  the  effects  of  the  variation  in 
the  dielectric  properties  of  the  medium  of  interest,  of  wavelength,  of 
the  angle  of  incidence,  and  of  transmit-recei ve  polarizations.  Back- 
scatter  theory  is  very  useful  as  an  aid  to  give  insight  into  the  under¬ 
standing  of  scattering  mechanisms. 

The  validity  of  the  theory  depends  very  heavily  on  the  mathematical 
model  that  is  used  to  describe  target  of  interest,  as  well  as  on  the 
approximations  that  may  be  necessary  to  be  made  to  obtain  some  solution 
to  the  problem.  Since  even  the  simplest  ground  surfaces  are  extremely 
difficult  to  describe  accurately,  a  true  mathematical  description  of 
the  ground  surface  or  target  volume  appears  impossible.  Therefore,  it 
is  necessary  to  measure  the  radar  return  from  natural  and  complex 
targets  and  then  to  use  theory  in  a  complementary  role;  i.e.,  as  an 
aid  in  interpreting  the  results  of  the  measurements,  and  in  suggesting 
how  they  may  be  extrapolated. 

5.1  Parashar's  Backscatter  Model 

Parashar  (1974)  made  an  attempt  to  describe  theoretically  the 
radar  backscatter  cross-section  for  sea  ice  by  considering  the  surface 
roughness  of  the  ice  and  the  amount  of  brine  entrapped  in  the  ice. 

Since  the  complex  permittivity  of  sea  ice  is  determined  by  both  tne 
properties  of  pure  ice  and  the  properties  and  amount  of  impurities 
present  in  the  ice,  Parashar  used  these  in  the  description  of  the 
vertical  profile  of  the  complex  permittivity.  Random  horizontal 
variations  of  the  complex  permittivity  were  also  taken  into  account. 

The  complex  permittivity  was  considered  to  have  a  prescribed  mean 
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vertical  profile  with  random  horizontal  variation  around  this  mean.  The 
relative  random  variation  was  assumed  to  be  small. 

The  wave  scattering  was  considered  to  have  components  both  ^rom  an 
irregular  air/ice  interface  and  from  a  continuously  varyi.g  inhomo¬ 
geneous  medium  (sea  ice).  The  salinity  profile  of  the  ice  was  approxi¬ 
mated  by  two  intersecting  line  segments  meeting  at  their  end  points.  It 
was  also  convenient  to  consider  the  ice  to  consist  of  two  regions,  joined 
at  the  breakpoint  in  the  salinity  profile.  Maxwell's  equations  were 
solved  using  the  small  perturbation  techniques  (Fung,  1969),  WKB  methods 
(Wait,  1962;  and  Brekhovski kh ,  1960),  and  variations  of  parameters 
(Hildebrand,  1962).  Calculations  were  performed  for  six  categories  of 
ice  to  obtain  radar  cross-sections  versus  incidence  angle  at  13.3  GHz 
and  400  MHz  for  VV  and  HH  polarizations.  Correlation  lengths  and 
variances  were  assumed  for  surface  roughness  and  for  the  horizontal 
inhomogeneity  within  the  ice  for  the  different  ice  categories  such  that 
good  agreement  existed  between  the  theoretical  and  experimental  results. 
Since  surface  truth  information  was  not  available,  parameters  such  as 
salinity,  temperature,  surface  roughness  and  complex  permittivity  were 
assumed  for  each  category  of  ice.  Experimental  and  theoretical  results 
showed  similar  relative  trends  from  category  to  category. 

5.2  Comparison  of  the  Results  of  the  Kansas  Experiments  with  Parashar's 

Theoretical  Model  and  the  Results  of  Others 

The  backscatter  model  correctly  predicted  that  the  returns  from 
multiyear  ice  should  be  greater  than  those  of  thick  first-year  ice 
(Figure  5-1).  The  model  indicated  a  difference  in  absolute  level  would 
be  on  the  order  of  2  dB  while  Parashar's  experimental  results  indicated 
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FIGURE  5-1 


A  comparison  of  Parashar's  Theoretical  Results  and  the 
Results  of  the  Recent  Kansas  Experiments  at  13  GHz,  VV. 


4  dB,  the  recent  Kansas  results  indicated  7  dB,  Gray's  results  indicated 
4  to  13  dB,  and  Rouse's  results  indicated  6  dB  (see  Figure  5-2).  The 
model  angular  responses  from  10°  to  60°  may  be  fit  with  a  straight  line. 
This  was  also  found  to  be  true  in  all  the  experimental  data.  The  model 
showed  returns  at  VV  polarization  were  consistently  higher  than  returns 
at  HH  polarization  and  their  slopes  were  significantly  different.  Recent 
Kansas  data  had  VV  returns  consistently  higher  than  HH  returns,  but  with 
similar  slopes.  The  model  indicated  that  the  response  of  multiyear  ice 
and  thick  first-year  ice  were  effectively  the  same  from  10°  to  60°  and 
only  differed  in  absolute  level.  Parashar's  and  Rouse's  data  indicated 
the  same,  but  our  recent  data  and  Gray’s  data  indicated  that  differences 
existed  in  the  slopes  on  the  two  responses.  The  absolute  levels  of  our 
recent  results  bounded  the  results  of  the  theoretical  model  and  the  abso¬ 
lute  levels  of  the  other  experimental  data  varied  significantly  from 
one  another  and  from  the  theoretical  model. 

Absolute  levels  and  the  angular  response  of  experimentally  measured 
radar  cross-sections  of  a  target  may  vary  from  sensor  to  sensor  as  was 
demonstrated  in  the  results  presented  above.  Inaccuracies  may  arise  due 
to  two  problem  areas.  First,  the  illumination  pattern  of  the  antennas 
must  be  known  in  great  detail  and  second,  extreme  care  must  be  taken  to 
calibrate  the  sensor  and  this  calibration  must  be  done  periodically. 

For  a  further  discussion  of  these  problem  areas  see  Chapter  1. 

5.3  An  Empirical  Model  of  Radar  Backscatter  from  Thick  First- 


The  angular  responses  of  the  radar  cross-sections  at  8  -  18  GHz 
frequencies  of  thick  first-year  sea  ice,  multiyear  sea  ice,  and  a  frozen 
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2  Comparison  of  multiyear  and  thick  first-year  ice 
scattering  cross-sections  (Experimental)  according 
to  Rouse  (1968),  Parashar,  et  al.  (1974),  and  Gray, 
et  al.  (1977). 


fresh  water  lake  with  an  underlying  layer  of  water  and  light  snow  cover 
were  found  to  be  reasonably  well  fit  with  a  straight  line  at  angles  of 
incidences  ranging  from  10°  to  70°.  The  straight  lines  were  obtained 
by  linear  regression  using  the  average  (a  magnitude  average  but  expressed 
in  dB)  scatter  obtained  from  the  multiple  looks  measured  at  a  given  site. 
The  scatter  was  found  to  fit  a  straight  line  for  both  like-  and  cross¬ 
polarizations.  The  average  correlation  coefficient  obtained  from  the 
regressed  data  for  HH  polarization  was  on  the  order  of  .97.  Using  the 
critical  absolute  value  test  to  indicate  the  significance  of  the  correla¬ 
tion  coefficient  results  in  a  critical  value  of  .87  for  a  1*  level  of 
significance.  By  comparing  the  correlation  coefficient  and  the  critical 
value  it  is  seen  that  the  null  hypothesis  that  there  is  zero  correlation 
in  these  linear  trends  is  not  justified.  An  empirical  model  has  been 
developed  to  describe  the  backscatter  for  these  three  types  of  ice  for 
the  8-18  GHz  range  of  frequencies,  angles  of  incidence  ranging  from  10° 
to  70°,  and  HH  polarization.  The  slopes  of  VV  and  HH  polarizations 
were  found  to  be  very  much  the  same  (both  visually  and  statistically) 
and  thus  a  model  for  VV  polarization  would  be  found  to  be  basically  the 
same. 

The  basic  regression  of  the  angular  response  of  the  data  is 
described  by: 

0°B  =  Y  +  Me  10°  <  e  <_  70°  (5-1) 

where: 

cr^g  is  the  differential  scattering  coefficient  expressed  in  dB 

Y  is  the  intercept  value 
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M  is  the  slope  in  dB/° 

and  a  is  the  angle  of  incidence  expressed  in  degrees. 

This  is  equivalent  to,  for  real  values  of  the  scattering  coefficient, 

-e/e„ 

ct°  =  aQe  0  10°  £  e  £  70°  (5-2) 

where: 

2  2 

a°  is  the  differential  scattering  coefficient  in  m  /m 
aQ  is  the  magnitude  constant 

and  9Q  is  the  e-folding  angle. 

For  our  observation  it  was  found  that  the  regression  must  take 
into  account  both  variations  in  return  due  to  angle  and  variations  due 
to  the  operating  frequency.  The  slopes  and  intercepts  obtained  from  the 
angular  responses  were  regressed  and  fit  with  a  straight  line.  These 
regression  equations  were 

M  =  A  +  Bf  8  £  f  <  18  GHz  (5-3) 

Y  =  D  +  Cf  8  £  f  £  18  GHz  (5-4) 

where 

A  is  the  intercept  value  of  the  frequency  response  of  the  slopes 
obtained  from  the  angular  responses 
8  is  the  slope  value  of  the  frequency  response  of  the  slopes 
obtained  from  the  angular  responses 
C  is  the  slope  value  of  the  frequency  response  of  the  intercepts 
obtained  from  the  angular  responses 
0  is  the  intercept  value  of  the  frequency  response  of  the  intercepts 
obtained  from  the  angular  responses 


and  f  is  the  value  of  the  frequency  ranging  from  8  to  18  (GHz). 

The  overall  expression  of  the  scattering  cross-section  is  then  described 
by  the  equation: 

8  <_  f  <_  18  GHz 

a°(dB)  =  (A+Bf )0  +  Cf  +  D  (5-5) 

10°  <  8  <  70° 

The  coefficients  for  the  three  ice  types  that  were  modeled  have  been 
included  in  Table  5-1.  Also  included  is  the  standard  error  associated 
with  the  model  when  compared  to  the  averaged  data.  This  error  was  2  dB 
or  less  and  was  considered  to  be  reasonable  because  it  maintained  the 
ability  to  discriminate  among  ice  types  and  because  the  error  was  less 
than  the  error  associated  with  the  data  and  the  averaged  data.  The 
standard  error  associated  with  the  data  was  around  3  dB  for  the  worst 
case.  For  a  comparison  of  the  three  ice  types  using  the  empirical  model 
and  a  display  of  the  standard  deviation  about  the  angular  response,  see 
Figure  5-3. 

A  statistical  test,  t-test,  was  performed  to  determine  if  the 
hypothesis  that  multiyear  sea  ice,  thick  first-year  sea  ice,  and  lake 
ice  were  statistically  different  was  true.  A  large  angle  of  incidence, 
60°,  was  chosen  for  this  test  because  the  separation  between  ice  types 
was  less,  thus  a  worst  case.  Data  at  9  GHz  (VV)  was  chosen  since  a 
sensor  with  these  parameters  is  very  much  of  interest.  It  was  found  that 
there  was  a  100%  confidence  that  multiyear  and  thick  first-year  sea  ice 
were  statistically  different.  This  was  also  found  to  be  true  for  multi¬ 
year  sea  ice  and  lake  ice.  There  was  a  94.7%  confidence  that  thick 
first-year  sea  ice  and  lake  ice  were  statistica1 ly  different. 


Regression  Model  Coefficient  for  Thick  First-Year  Sea  Ice,  Multiyear  Sea  Ice 
and  Fresh  Water  Lake  Ice  with  an  Underlying  Layer  of  Water  and  Little  Snow  Cover  for 


5.4  A  Qualitative  Description  of  Backscatter 


A  qualitative  description  of  the  backscatter  properties  of  thick  first- 
year  sea  ice,  multiyear  sea  ice,  and  fresh  water  lake  ice  at  Ku-X-Band 
frequencies  has  been  made.  From  this  description  the  existing  backscatter 
theory  may  be  modified  and  extended  where  needed.  The  scattering  medium 
will  be  considered  to  consist  of  three  major  regions.  The  first  is  a 
homogeneous  layer  of  air;  the  second  is  a  layer  of  ice  which  will  be 
described  according  to  its  surface  roughness,  its  inhomogenuities,  and 
its  dielectric  properties;  and  a  homogeneous  layer  of  sea  water. 

A  family  of  radar  cross-section  curves,  Figure  5-3,  has  been  generated 
to  illustrate,  with  many  qualifications,  the  complex  nature  of  the  inter¬ 
action  of  surface  roughness  and  the  volume  scatter  properties  of  the  ice 
medium  on  angular  response.  Five  conditions  were  considered.  They  are 
backscatter  due  to  volume  scatter  alone,  to  a  dominant  volume  scatter  with 
a  contribution  of  scatter  from  a  slightly  rough  (<x)  surface,  to  scatter 
from  a  slightly  rough  surface  with  some  volume  scatter,  to  scatter  from 
a  slightly  rough  surface,  and  to  scatter  from  a  very  slightly  rough 
surface.  From  this  set  of  curves  speculation  into  the  backscatter  properties 
of  ice  will  be  made.  In  Figure  4.5-33,  the  angular  response  of  three 
categories  of  ice  is  illustrated.  In  our  discussion,  backscatter  will  be 
described  in  terms  of  the  contribution  of  volume  scatter,  surface  scatter, 
or  the  interaction  of  the  two. 

First  we  consider  multiyear  ice.  Backscatter  is  speculated  to 
consist  of  a  significant  component  of  volume  scatter  and  a  component 
of  surface  scatter.  Volume  scatter  is  suggested  due  to  the  gradually 
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FIGURE  5-3 


Scattering  coefficients  of  multiyear,  thick  first-year 
and  lake  ice  at  9  GHz,  HH,  using  an  empirical  model 
developed  from  the  recent  Kansas  experiments.  Dashed 
lines  indicate  +  s...... 


A  family  of  radar  cross-sections  which  illustrate  the 
interaction  of  surface  roughness  and  the  volume  scatter 
properties  of  ice. 
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decreasing  angular  response.  Returns  at  large  angles  fall  off  too 
rapidly  to  be  volume  scatter  alone.  The  addition  of  a  rough  surface 
has  the  effect  of  decreasing  the  isotropic  effect  of  volume  scatter 
(associated  with  a  smooth  surface)  which  results  in  a  more  rapidly 
decreasing  angular  response.  The  speculation  that  volume  scatter  is 
a  significant  component  is  seen  to  be  justifiable  when  the  attenuation 
coefficient  and  skin  depth  are  examined.  According  to  the  literature 
[60,61],  the  attenuation  coefficient  is  in  the  range  of  10  dB  or  less 
for  1"  salinity  ice,  -5°  C,  and  10  GHz.  The  depth  at  which  the  electric 
field  in  the  ice  has  been  attenuated  to  37%  of  its  original  value  is 
called  the  skin  depth  and  is  1.2  meter  for  the  conditions  given  above. 

Two  skin  depths  shall  be  considered  to  be  a  reasonable  choice  from 
which  to  consider  the  extent  of  the  scattering  volume  into  the  ice. 

The  power  at  this  point  is  2%  of  the  power  transmitted  into  the  ice. 

Therefore,  the  scattering  volume  for  the  above  conditions  would  extend 
into  tne  multi year  ice  2.4  m.  Surface  truth  shows  dense  distributions 
of  air  bubbles  of  a  size  less  than  a  wavelength.  These  serve  as  the 
scattering  sources  in  the  volume  and  indicate  the  capability  of  the 
ice  to  have  a  large  component  of  volume  scatter. 

It  is  speculated  that  the  radar  cross-section  of  thick  first-year 
ice  is  dominated  by  surface  scatter  with  some  contribution  of  volume 
scatter.  The  rapidly  decaying  angular  response  near  nadir  angles  is 
a  shape  that  is  characteristic  of  scatter  from  a  slightly  rough  surface. 

However,  a  hump  occurs  at  the  middle  angles  and  its  existence  is 
explained  by  including  a  component  of  volume  scatter.  The  above 
speculation  about  backscatter  from  this  type  of  ice  is  justified  by 
noting  that  the  surface  has  a  roughness  that  is  on  the  order  of  less 

i 
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than  a  wavelength,  the  attenuation  coefficient  for  8%  salinity  ice  at 
-5°  C  and  10  GHz  is  300-500  dB/m,  and  the  skin  depth  is  2-3  cm  [60,62]. 
Therefore,  the  scattering  volume  is  small  (the  top  4-6  cm  of  the  ice 
sheet)  and,  as  seen  from  the  surface  truth  photos,  bubbles  and  inclusions 
on  the  order  of  much  less  than  a  wavelength  exist  in  this  upper  layer 
of  the  ice  sheet  and  serve  as  scattering  centers. 

Lake  ice  has  an  even  lower  return  than  thick  first-year  ice. 
Speculation  is  that  the  surface  scatter  contribution  is  on  the  same 
order  of  magnitude  as  that  of  the  volume  scattering  component.  This 
speculation  considers  the  very  uniform  and  decreasing  angular  response. 
The  surface  is  much  smoother  than  the  surface  of  thick  first-year  ice 
and  roughness  is  on  an  order  of  much-less-than-a-wavelength.  Therefore, 
the  absolute  contribution  due  to  surface  scatter  will  be  less.  The 
attenuation  coefficient  is  on  the  order  of  1  dB/m  and  a  skin  depth  of 
8.7  m.  Therefore,  the  scattering  volume  extends  to  the  ice/water  inter¬ 
face.  Small  air  bubbles,  much  less  than  a  wavelength  in  size,  are  found 
throughout  the  ice.  Long  drain  tubes  with  diameters  on  the  order  of 
.5  mm  and  lengths  on  the  order  of  wavelengths  are  found  near  the  bottom. 
However,  in  comparison  with  multiyear  and  thick  first-year  ice  the  popu¬ 
lation  of  scattering  centers  is  light.  The  reflection  from  the  ice/water 
interface  also  serves  to  enhance  the  volume  scattering  effect. 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Results 

The  ability  to  use  radar  to  discriminate  ice  types  has  further 
been  demonstrated.  Thick  first-year,  multiyear,  and  pressure  ridged 
sea  ice  and  fresh  water  lake  ice  have  been  investigated  at  1.5  and 
from  8-18  GHz,  at  angles  of  incidence  from  10°  to  70°  (with  respect  to 
nadir),  and  at  antenna  transmi t-recei ve  polarizations  of  VV  and  VH  for 
1.5  GHz  and  VV,  VH,  HH,  and  HV  for  8-18  GHz. 

The  L-Band  radar  was  found  to  have  the  ability  to  discriminate 
among  fresh  water  ice,  the  thick  categories  of  sea  ice,  and  pressure 
ridged  sea  ice.  At  angles  greater  than  60°  discrimination  between 
thick  first-year  and  multiyear  ice  may  be  possible,  especially  at 
cross-polarization.  What  was  found  in  these  experiments  confirms  what 
had  been  concluded  from  earlier  SAR  imagery.  The  practical  use  of  L- 
Band  radar  would  be  in  topographical  mapping. 

The  Ku-X-Band  radar  was  able  to  discriminate  among  the  ice  types 
investigated.  The  lowest  frequency,  9  GHz,  was  found  to  provide  the 
greatest  separation  between  two  ice  types  of  major  concern,  thick 
first-year  and  multiyear  sea  ice.  Cross-polarization  provided  a  slightly 
greater  separation  then  VV  and  HH.  Significant  levels  of  separation 
existed  between  angles  of  incidence  from  15°  to  70°.  Pressure  ridges 
are  also  of  great  concern.  Since  they  are  linear  features  of  high 
return  on  a  thick  fi rst-year  background,  they  should  be  recognizable 
due  to  these  characteristics  in  spite  of  possible  confusion  and  an 
absolute  scattering  level  basis  with  multiyear  ice. 
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Ku-X-Band  angular  responses  from  10°  to  70°  were  regressed  with 
a  straight  line.  An  average  correlation  coefficient  of  .97  (for  1  % 
confidence  the  correlation  coefficient  for  our  measurements  would 
need  to  be  greater  than  .89)  shows  the  goodness  of  the  fit.  The  inter¬ 
cepts  and  slopes  from  this  straight  line  were  also  regressed  with 
straight  lines  and  the  results  were  used  in  an  empirical  model  which 
describes  the  angular  and  frequency  variations  for  thick  first-year 
and  multiyear  sea  ice  and  lake  ice.  The  standard  deviation  between 
the  model  and  averaged  scattering  cross-sections  was  at  worst  2  dB. 

The  individual  measurements  had  a  worst  standard  deviation  of  3  dB 
about  the  average  scattering  cross-sections.  Thus  the  model  effectively 
describes  the  data. 

6.2  Recommendations 

There  is  a  great  need  to  investigate  further  the  radar  backscatter 
properties  of  sea  ice.  The  experiments  up  to  this  point  have  surveyed, 
in  a  sense,  what  may  be  done  with  radar.  The  1977  and  1978  Kansas 
experiments  were  the  first  of  what  hopefully  will  be  many  carefully 
staged  attempts  to  describe  quantitatively  the  backscatter  properties 
of  sea  ice.  In  these  experiments  an  abundance  of  surface  truth  infor¬ 
mation  has  been  collected  which  describes  the  ice  under  investigation. 
This  information  will  be  used  to  develop  an  ice  model  from  which  a 
second  generation  backscatter  model  may  evolve.  Further  work  desperately 
needs  to  be  done  to  pinpoint  the  locations  and  elements  in  the  sea  ice 
medium  which  contribute  to  the  scatter.  Recording  of  the  intermediate- 
frequency  returns  from  the  target  and  processing  these  returns  using 
Fast-Fourier-transform  techniques  have  been  proposed  to  begin  this  task. 
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Study  of  the  effects  of  temperature  on  backscatter  return  is  also 
needed.  This  must  include  measurement  programs  in  the  growing  and 
melting  seasons.  These  tend  to  be  inherently  difficult  logistically, 
but  are  necessary  in  a  complete  description  of  the  radar  return  of 
sea  ice. 

The  1978  experiment  was  the  first  experiment  to  consider  quanti¬ 
tatively  the  role  of  the  snow  layer  as  a  contributor  to  backscatter. 
Differences  between  returns  from  ice  with  a  snow  cover  and  bare  ice 
were  in  the  3  dB  range  for  thick  first-year  ice  and  up  to  8  dB  for  lake 
ice.  These  effects  need  to  be  considered  in  the  study  of  the  signa¬ 
ture  of  multiyear  ice  where  the  ice  is  almost  pure  and  snow  depths 
vary  from  none  to  many  tens  of  cm. 

The  9  GHz  frequency  was  found  to  be  the  best  for  discrimination. 
But  what  about  2-8  GHz?  It  follows  that  these  frequencies  should  also 
be  considered  in  the  puzzle  since  the  lower  X-Band  frequencies  were 
found  to  work  so  well. 

Comparison  of  the  measurements  taken  with  the  surface-based  radar, 
the  airborne  SLAR  and  the  SARs  would  be  useful.  Besides  serving  to 
give  a  calibration  to  the  airborne  imagery  it  would  indicate  the 
peculiarities  that  may  arise  due  to  the  differences  in  what  each  radar 
sees. 

What  has  yet  to  be  done  and  would  be  extremely  valuable  is  the  use 
of  airborne  imaging  radars  in  conjunction  with  a  highly  mobile  scatter- 
ometer  and  ground  truth  team.  The  scheme  of  this  type  of  complementary 
experiment  is  the  use  of  an  airborne  radar  to  provide  real  time  images, 
pinpoint  areas  of  greatest  interest  from  these  images,  fly  the  surface 
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truth  team  to  the  sites  of  interest,  investigate  the  physical  proper¬ 
ties  of  these  sites,  and  have  a  mobile  calibrated  radar  make  backscatter 
measurements.  A  fundamental  problem  has  existed  in  past  experiments 
in  that  imagery  has  been  processed  after  the  experiment  is  over.  Be¬ 
cause  of  this,  sites  that  prove  most  puzzling  and  may  contain  an  abun¬ 
dance  of  information  remain  unseen  by  surface  parties. 

The  highly  mobile  calibrated  radar  that  has  been  referred  to 
above  has  been  developed  by  the  University  of  Kansas  (HELOSCAT).  It 
is  mountable  to  a  helicopter  and  at  present  operates  at  frequencies  of 
8-18  GHz,  HH  polarization,  and  incidence  angles  of  20°,  40°,  and  60°. 


6.3  Sensor 

Description  of  an  ice  reconnaissance  radar  follows  from  the  analy¬ 
sis  of  the  data  that  is  presently  available.  The  operating  frequency 
would  be  at  the  lowest  X-Band  frequencies  (example:  9  GHz),  operate  at 
VV  polarization,  and  at  angles  of  incidence  greater  than  25°. 

The  lowest  X-Band  frequencies  were  chosen  due  to  the  frequency 
responses  measured  in  the  Kansas  experiments.  Indications  are  that 
separation  among  sea  ice  types  is  greater  at  this  lower  X-Ku-Band  fre¬ 
quency.  VV-polarization  was  selected  according  to  Parashar's  theoretical 
backscatter  model  and  the  results  of  the  recent  Kansas  experiments. 

VV  returns  allow  the  same  degree  of  discrimination  as  HH  but  radar 
cross-sections  are  slightly  higher  (recent  experiments)  or  significantly 
higher  (Parashar)  in  absolute  level.  Cross-polarization  would  allow 
a  couple  of  dB  of  added  separation  to  a  range  of  separation  that  is 
already  large  for  thick  first-year  and  multiyear  ice,  but  this  is  at 


the  expense  of  working  with  cross-sections  that  are  on  the  order  of 
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10  dB  lower.  Nothing  can  be  said  at  this  time  to  indicate  that  a 
real  benefit  exists  in  operating  at  cross-pol arization;  however, 

Parashar  and  CCRC  data  indicate  a  possible  advantage  for  thin  ice  types 
not  studied  in  the  recent  experiments.  Angles  of  incidence  greater 
than  25°  were  chosen  by  considering  both  recent  results  and  Parashar' s 
theoretical  and  experimental  results.  Least  confusion  in  discrimination 
occurs  at  angles  away  from  the  region  where  returns  from  different  cate¬ 
gories  of  ice  intersect  and  change  trends  (15°  to  20°).  Angles  between 
nadir  and  this  region  of  overlap  also  may  permit  discrimination.  How¬ 
ever,  open  water  has  a  high  return  at  these  angles  and  the  effects  of 
wind  speeds  on  its  backscatter  return  must  be  considered  in  the  problem. 
At  angles  greater  than  25°  the  angular  responses  more  or  less  flatten 
out.  The  groups  of  categories  of  ice  that  would  be  discriminable 
according  to  the  results  at  hand  are  multiyear  (>2m),  thick  first- 
year  (l-2m),  new  and  thin  ice  (0+-18cm),  and  thick  young  and  thin 
ice  (18-100cm).  Open  water  in  many  instances  would  produce  very  low 
returns  at  these  angles  and  be  discriminable. 
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